ON THE BOLTZMANN-GRAD LIMIT 
FOR THE TWO DIMENSIONAL PERIODIC LORENTZ GAS 



EMANUELE CAGLIOTI AND FRANQOIS GOLSE 

Abstract. The two-dimensional, periodic Lorentz gas, is the dynamical sys- 
tem corresponding with the free motion of a point particle in a planar system 
of fixed circular obstacles centered at the vertices of a square lattice in the Eu- 
clidian plane. Assuming clastic collisions between the particle and the obsta- 
cles, this dynamical system is studied in the Boltzmann-Grad limit, assuming 
that the obstacle radius r and the reciprocal mean free path are asymptoti- 
cally equivalent small quantities, and that the particle's distribution function 
is slowly varying in the space variable. In this limit, the periodic Lorentz 
gas cannot be described by a linear Boltzmann equation (see [F. Golse, Ann. 
Fac. Sci. Toulouse 17 (2008), 735-749]), but involves an integro-differential 
equation conjectured in [E. Caglioti, F. Golse, C.R. Acad. Sci. Ser. I Math. 
346 (2008) 477-482] and proved in [J. Marklof, A. Strombergsson, preprint 
arXiv:0801.0612], set on a phase-space larger than the usual single-particle 
phase-space. The main purpose of the present paper is to study the dynami- 
cal properties of this integro-differential equation: identifying its equilibrium 
states, proving a H Theorem and discussing the speed of approach to equilib- 
rium in the long time limit. In the first part of the paper, we derive the explicit 
formula for a transition probability appearing in that equation following the 
method sketched in [E. Caglioti, F. Golse, loc. cit.]. 



1. The Lorentz gas 

The Lorentz gas is the dynamical system corresponding with the free motion of a 
single point particle in a system of fixed spherical obstacles, assuming that collisions 
between the particle and any of the obstacles are elastic. This simple mechanical 
model was proposed in 1905 by H.A. Lorentz [17] to describe the motion of electrons 
in a metal — see also the work of P. Drude [10] 

Henceforth, we assume that the space dimension is 2 and restrict our attention 
to the case of a periodic system of obstacles. Specifically, the obstacles are disks 
of radius r centered at each point of Z 2 . Hence the domain left free for particle 
motion is 

(1.1) Z r = {x e R 2 | dist(x, Z 2 ) > r} , where < r < \. 

Throughout this paper, we assume that the particle moves at speed 1. Its tra- 
jectory starting from x G Z r with velocity u> G S 1 at time t = is denoted by 
1 1 ^ (X r ,n r )(t;x,u}) € R 2 x S 1 . One has 

(1.2) | 2 o rW ' whenever X r (t) e Z r , 
while 

^ {Zllo) ZntmAt-0), whenever 
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In the system above, we denote ' = ^, and lZ[X r (t)} is the specular reflection on 
dZ r at the point X r (t) = X r (t ± 0). 

Next we introduce the Boltzmann-Grad limit. This limit assumes that r ^ 1 
and that the initial position x and direction uj of the particle are jointly distributed 
in Z r x S 1 under some density of the form f m (rx,td) — i.e. slowly varying in x. 
Given this initial data, we define 

(1.4) f r (t,x,uj) :— f m (rX r (—t/r;x,uj) 7 £l r (—t/r;x 7 uj)) whenever x G Z r . 

In this paper, we are concerned with the limit of f r as r — > + in some sense to be 
explained below. In the 2-dimensional setting considered here, this is precisely the 
Boltzmann-Grad limit. 

In the case of a random (Poisson), instead of periodic, configuration of obstacles, 
Gallavotti [11] proved that the expectation of f r converges to the solution of the 
Lorentz kinetic equation 



(1.5) 



(d t +w-'V x )f(t,x,u>) = / (f(t,x,u)-2(u> ■ n)n)-f(t,x,uj))(u) ■ n) + dn, 
Js 1 

f 



for all t > and (x,u) e R 2 x S 1 . Gallavotti's remarkable result was later gener- 
alized and improved in [23, 3]. 

In the case of a periodic distribution of obstacles, the Boltzmann-Grad limit of 
the Lorentz gas cannot be described by the Lorentz kinetic equation (1.5). Nor can 
it be described by any linear Boltzmann equation with regular scattering kernel: 
sec [12, 14] for a proof of this fact, based on estimates on the distribution of free 
path lengths to be found in [4] and [15]. 

In a recent note [7], we have proposed a kinetic equation for the limit of f r as 
r — > + . The striking new feature in our theory for the Boltzmann-Grad limit of 
the periodic Lorentz gas is an extended single-particle phase-space (see also [13]) 
where the limiting equation is posed. 

Shortly after our announcements [13, 7], J. Marklof and A. Strombergsson inde- 
pendently arrived at the same limiting equation for the Boltzmann-Grad limit of 
the periodic Lorentz gas as in [7]. Their contribution [19] provides a complete rig- 
orous derivation of that equation (thereby confirming an hypothesis left unverified 
in [7]), as well as an extension of that result to the case of any space dimension 
higher than 2. 

The present paper provides first a complete proof of the main result in our note 
[7] . In fact the method sketched in our announcement [7] is different from the one 
used in [20] , and could perhaps be useful for future investigations on the periodic 
Lorentz gas in 2 space dimensions. 

Moreover, we establish some fundamental qualitative features of the equation 
governing the Boltzmann-Grad limit of the periodic Lorentz gas in 2 space dimen- 
sions — including an analogue of the classical Boltzmann H Theorem, a description 
of the equilibrium states, and of the long time limit for that limit equation. 

We have split the presentation of our main results in the two following sec- 
tions. Section 2 introduces our kinetic theory in an extended phase space for the 
Boltzmann-Grad limit of the periodic Lorentz gas in space dimension 2. Section 3 
is devoted to the fundamental dynamical properties of the integro-differential equa- 
tion describing this Boltzmann-Grad limit — specifically, we present an analogue 
of Boltzmann's H Theorem, describe the class of equilibrium distribution functions, 
and investigate the long time limit of the distribution functions in extended phase 
space that are solutions of that integro-differential equation. 
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2. Main Results I: The Boltzmann-Grad Limit 

Let (x, uj) G Z r x S 1 , and define < to < t\ < . . . to be the sequence of collision 
times on the billiard trajectory in Z r starting from x with velocity uj. In other 
words, 

(2.1) {tj | j € N} = {t £ H* + | X r (t; x, uj) G dZ r } . 
Define further 

(2.2) (xj,u} j ):=(X r (t j +0;x,uj),n r (t j +0;x,uj)), j>0. 

Denote by n x the inward unit normal to Z r at the point x £ dZ r , and consider 
L± = {(x,lj) e dZ r x S 1 | ±u -n x > 0}, 
U '" >) T° r = {(x,u) edZ r xS 1 \uj-n x = 0}. 

Obviously, (xj,u)j) G L+ U T° for each j > 0. 

2.1. The Transfer Map. As a first step in finding the Boltzmann-Grad limit of 
the periodic Lorentz gas, we seek a mapping from L+ U Tj? to itself whose iterates 
transform (xq,uo) m to the sequence (xj,L)j) defined in (2.2). 

For (x,lo) e (Z r x S 1 ) U L+ U let T r (x,w) be the exit time defined as 

(2.4) r r (x, uj) = w£{t >0\x + tuje dZ r } . 

Also, for (x,uj) G U Tjl, define the impact parameter h r {x,uj) as on Figure 1 by 

(2.5) h r (x, u> ) = sin^', n x ) . 

Denote by (T+ UT°)/Z 2 the quotient of T+ U Tj! under the action of Z 2 by 
translations on the x variable. Obviously, the map 

(2.6) (r+ U T°)/Z 2 3 (x, uj) ^ (h r (x, uj),uj) G [—1,1] x s 1 

coordinatizes (r+ U r°)/Z 2 , and we henceforth denote by Y r its inverse. For r G 
]0, i[, we define the transfer map (see Figure 2) 

T r : [-1,1] x S 1 ->• R; x [-1,1] 

by 

(2.7) T r {h',u) = (2r Tr {Y r {h' 1 uj)) 1 h r ((X ri n r ){T r {Y r (h' 1 uj))±0;Y r {h' 1 uj)))). 
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Figure 2. The transfer map 



Up to translations by a vector of Z 2 , the transfer map T r is essentially the sought 
transformation, since one has 

(2.8) T r (h r (xj,ujj),ujj) = (2rr r (xj,ujj),h r (xj + i,ujj)) , for each j > 0, 
and 



-i = R[tt — 2 arcsin(/i r (:Ej + i , , for each j > . 



(2.9) 
The notation 

(2.10) R[9] designates the rotation of an angle 9. 
Notice that, by definition, 

h r {x j+ i,ujj) = h r (x j+ i,ujj + i) . 

The theorem below giving the limiting behavior of the map T r as r — > + was 
announced in [7]. 

Theorem 2.1. For each $ e C C (R* + x [-1, 1]) and each h! € [-1, 1] 

(2.11) - — r / $(TJti,u)) >• / / $(S 7 h)P(S,h\ti)dSdh 

|lne|J e r J J_ x 

a.e. in uj £ S 1 as e — > + , where the transition probability P(S, h\h')dSdh is given 

by the formula 

(2.12) 

P(S,h\h>) = ^((Sv) A (l-S)+ + (riS-\l-S\) + 



with the notation 
and 



+ ((s a (i + \Sr,) (is + isc) v i) + 

+ ((S - iSry) A 1 - (|S + \SQ V (1 - iSr?)) + ) , 
C = i|/i + /i'|, ri = \\h-ti\, 



a Ab — inf (a, b) , a V & = sup(a, 6) . 
Equivalently, for each (S,h,h') € R+x] — l,l[x] — 1,1[ suc/i i/iai < ft, 



(2.13) 



P(S,ft|ft') = £ 1A 



1 



ft - ft' V S 



(1 + ft') 



The formula (2.12) implies the following properties of the function P. 
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Corollary 2.2 (Properties of the transition probability P(S, h\h').). The function 
(S, h, h') i y P(S, h\h') is piecewise continuous on R + x [—1, 1] x [—1, 1]. 

1) It satisfies the symmetries 
(2.14) 

P(S, h\ti) = P(S, h'\h) = P(S, -h\ - ti) for a.e. h, h! E [-1, 1] and S > , 
as well as the identities 

oo pi 

J P{S, h\h')dSdh = 1 , for each h' E [1, 1] , 

J P(S, h\h')dSdh' = 1 , for each h E [1, 1] . 

2) The transition probability P(S,h\h') satisfies the bounds 

(2.16) Q< P{Sl h\h')<^l l+w<i 

for a.e. h, h! g] — 1, 1[ such that \h'\ < h and all S > 4. Moreover, one has 

/l rl AO 
J P{S,h\h')dhdh' < -2^3 , S*>4. 

As we shall see below, the family T r (h',u>) is wildly oscillating in both h! and u> 
as r — > + , so that it is somewhat natural to expect that T r converges only in the 
weakest imaginable sense. 

The above result with the explicit formula (2.12) was announced in [7]. At the 
same time, V.A. Bykovskii and A.V. Ustinov 1 arrived independently at formula 
(2.13) in [5]. That formulas (2.12) and (2.13) are equivalent is proved in section 6.2 
below. 

The existence of the limit (2.11) for the periodic Lorentz gas in any space dimen- 
sion has been obtained by J. Marklof and A. Strombergsson in [18], by a method 
completely different from the one used in the work of V.A. Bykovskii and A.V. 
Ustinov or ours. However, at the time of this writing, their analysis does not seem 
to lead to an explicit formula for P(s, h\h') such as (2.12)-(2.13) in space dimension 
higher than 2. 

Notice that J. Marklof and A. Strombergsson as well as V.A. Bykovskii and A.V. 
Ustinov obtain the limit (2.11) in the weak-* L°° topology as regards the variable 
io, without the Cesaro average over r, whereas our result, being based on Birkhoff 's 
crgodic theorem, involves the Cesaro average in the obstacle radius, but leads to a 
pointwise limit a.e. in co. 

In space dimension 2, [20] extends the explicit formula (2.12)-(2.13) to the case 
of interactions more general than hard-sphere collisions given in terms of their 
scattering map. The explicit formula proposed by Marklof-Strombergsson in [20] 
for the transition probability follows from their formula (4.14) in [18], and was 
obtained independently from our result in [7]. 

Another object of potential interest when considering the Boltzmann-Grad limit 
for the 2-dimcnsional periodic Lorentz gas is the probability of transition on im- 
pact parameters corresponding with successive collisions, which is essentially the 
Boltzmann-Grad limit of the billiard map in the sense of Young measures. Ob- 
viously, the probability of observing an impact parameter in some infinitesimal 
interval dh around h for a particle whose previous collision occured with an impact 
parameter h' is 

/•OO 

U(h\h') = / P(S,h\h')dS. 
Jo 




We are grateful to J. Marklof for informing us of their work in August 2009. 
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Explicit formula for H(h\h') 



For \h'\ < h < 1, one has 



1 , 1 + h 



(2-18) ^ = ^h^ l «l + h > 

Besides, the transition probability U(h\h') satisfies the symmetries inherited from 
P(S,h\h'): 

(2.19) U(h\h') = TL(h'\h) = U(-h\ - ti) for a.e. h, h! e [-1, 1] . 

2.2. 3-obstacle configurations. Before analyzing the dynamics of the Lorentz 
gas in the Boltzmann-Grad limit, let us describe the key ideas used in our proof of 
Theorem 2.1. 

We begin with an observation which greatly reduces the complexity of billiard 
dynamics for a periodic system of obstacles centered at the vertices of the lattice 
Z 2 , in the small obstacle radius limit. This observation is another form of a famous 
statement about rotations of an irrational angle on the unit circle, known as "the 
three-length (or three-gap) theorem" , conjectured by Steinhaus and proved by V.T. 
Sos [22] — see also [24]. 

Assume weS 1 has components uj 1 , uj 2 independent over Q. Particle trajectories 
leaving, say, the surface of the obstacle centered at the origin in the direction uj will 
next collide with one of at most three, and generically three other obstacles. 

Lemma 2.3 (Blank-Krikorian [1], Caglioti-Golse [6]). Let < r < \, and uj e S 1 
be such that < uj 2 < uj 1 and uj 2 /uj 1 ^ Q. Then, there exists (q,p) and (q,p) in 
Z 2 such that 

0<q<q, qp — qp = a G {±1} 
satisfying the following property: 

{x + T r (x, lu)uj | \x\ = r and x ■ lu > 0} 

C dD((q,p), r) U dD((q,p), r) U dD{{q + q,p + p), r) , 
where D(xo,r) designates the disk of radius r centered at xq. 

The lemma above is one of the key argument in our analysis. 
To go further, we need a convenient set of parameters in order to handle all these 
3-obstacle configurations as the direction ui runs through S 1 . 
For uj as in Lemma 2.3, the sets 

{x + tco | |a;| = r , x ■ uj > , x + t v (x,uj)uj € dD((q,p),r) , t e R} 

and 

{x + toj | |a;| = r , x ■ uj > , x + r r (x, uj)uj S dD((q,p),r) , t e R} 

are closed strips, whose widths are denoted respectively by a and b. The following 
quantities are somewhat easier to handle: 

(2-20) Q = ^, Q = ^, A=±, B = ± 

uj l uj l 2r 2r 

(see Figure 3 for the geometric interpretation of A,B,Q and Q), and we shall 
henceforth denote them by 

(2.21) Q(uj,r) , Q(uj,r) , A(uj,r) , B{w,r) , together with a(co, r) 

whenever we need to keep track of the dependence of these quantities upon the 
direction uj and obstacle radius r — we recall that 

a(uj,r) =qp-pqe {±1} . 
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Figure 3. Example of a 3-obstacle configurations, and the pa- 
rameters A, B, Q and Q. Here e = 2r/w\. 



Lemma 2.4. Let < r < \, and uj € S 1 be such that < u 2 < uj 1 and lo 2 /u 1 ^ Q. 

Then, one has 

j < A(w, r) , B(u, r) , A(w, r) + B(w, r) < 1 , 
| < Q(u,r) < Q(fi,r), 



(2.22) 
and 



Q(u, r)(l - i4(w, r)) + Q(w, r)(l - S(w, r)) = 1 . 
This last equality entails the bound 

(2-23) 0< Q( OJj r)<— — 1 —- -<1. 

2 — A(lj, r) — B(w, r) 

Therefore, each possible 3-obstaclc configuration corresponding with the di- 
rection to and the obstacle radius r is completely determined by the parameters 
(A,B,Q,a)(u,r)e [0,l] 3 x{±l}. 

The proof of Theorem 2.1 is based on the two following ingredients. 

The first is an asymptotic, explicit formula for the transfer map T r in terms of 
the parameters A, B, Q, a defined above. 

Proposition 2.5. Let < r < \, and uj e S 1 be such that < uj 2 < uj 1 and 

uj 2 /uj 1 $l Q. Then, for each h' G [—1, 1] and each r e]0, \[, one has 

T r (ti,Uj) = TA(u,r),B(u,r),Q(u,r),<r(u,r)(h') + (O(r 2 ),0) , 

in the limit as r — > + . In the formula above, the map TA,B,Q,a * s defined for each 
(A,B,Q,a) £ [0,1] 3 x {±1} in the following manner: 

Ta,b,qAH) = (Q> h ' - 2ff (! - A )) if *K € [1 - 2A, 1] , 

(2.24) T A ,B,Q,a(ti) = {Q,h' + 2a{l-B)) if ah' G [-1,-1 + 2B] , 
T A ,B,C)A h ) = {Q + Q,h' + 2a(A - B)) otherwise. 

For ui = (cos 9, sin 9) with arbitrary 6 <G R, the map h' T r (h',uj) is computed 
using Proposition 2.5 by using the symmetries in the periodic configuration of 
obstacles as follows. Set 9 = 9 — ra-| with m = [§(#+§)] (where [z] is the integer 
part of z) and let Co = (cos 9, sin 9). Then 

(2.25) T r (h',u>) = (s,h) , where (s,sign(tan0)ft) = T r (sign(tan(9)/j', w) . 



8 



E. CAGLIOTI AND F. GOLSE 



The second ingredient in our proof of Theorem 2.1 is an explicit formula for the 
limit of the distribution ofwi4 (A, B, Q, a)(ui, r) as r — » + in the sense of Cesaro 
on the first octant S+ of S 1 . 

Proposition 2.6. Let F be any bounded and piecewise continuous function defined 
on the compact K = [0, l] 3 x {±1}. Then 



1 f L/i dr 
- — - / F(A(u,r),B(uj,r),Q(uj,r),a(uj,r)) — 

(2.26) l ln£ l^ r 

-> / F{A,B,Q,a)dfi(A,B,Q,a) 

JK 

a.e. inueS| as e — > + 7 w/iere (U is the probability measure on K given by 
(2.27) 

i / a t-> \ ^ i 1 dAdBdQ 
d^{A,B,Q 1 a)^—^\ 0<A< i\ 0<B< i-Ai-o < Q < ^-^- B y-A + <V=-i) • 

This result is perhaps more transparent when stated in terms of the new param- 
eters A 1 B' = j^j,Q, o instead of the original A, B 1 Q, a: an elementary change of 
variables in the integral on the right hand side of (2.26) shows that 



F(A,B,Q,o-)dfi(A,B,Q,a) 



f F(A,(l-A)B / ,Q,a)dn'(A,B',Q,a). 



where 



djj!{A, B',Q,a) 



= ■^■1o< j 4<i1o<B'<i1 <q< 1+(1 _a )(1 _ b i- ) dAdB 1 dQ{8 a= i + 5 a= -i) . 

In other words, the new parameters A, B' ,Q and a are uniformly distributed over 
the maximal domain compatible with the bounds (2.22) and (2.23). 

The first part of Theorem 2.1 follows from combining the two propositions above; 
in particular, for each h' € [—1, 1], the transition probability P(S, h\h')dSdh is ob- 
tained as the image of the probability measure fj, in (2.27) under the transformation 
(A,B,Q,a) i— > T AiBtQtir (ti). 



2.3. The Limiting Dynamics. With the parametrization of all 3-obstacle config- 
urations given above, we return to the problem of describing the Boltzmann-Grad 
limit of the Lorentz gas dynamics. 

Let (x,u) G Z r x S 1 , and let the sequence of collision times (tj)j>o, collision 
points (xj)j>o and post-collision velocities {u>j)j>o be defined as in (2.1) and (2.2). 
The particle trajectory starting from x in the direction w at time t = is obviously 
completely defined by these sequences. 

As suggested above, the sequences {tj)j>o, (xj)j>o and (wj)j>o can be recon- 
structed with the transfer map, as follows. 

Set 

t = T r (x,Cj) , 
In = X + T r (x,U))bJ , 
( 2 - 28 ) h h ( \ 

h = h r (x 0l uj) , 

uj a = R[n — 2 arcsin(/i )]w • 
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We then define the sequences (tj)j>o, (xj)j>o inductively, in the following manner: 
(2s j+1 ,h j+1 ) = T r (hj,Uj) , 
1 

tj+i — tj + -Sj+i , 

(2.29) \ 

Xj+l = Xj + -Sj+lUj , 

= R[n — 2arcsin(/i J+ i)]a;j , 

for each j > 0. 

If the sequence of 3-obstacle configuration parameters 

= ((A,B,Q,a)(w j ,r)) j > 

converges (in some sense to be explained below) as r — > + to a sequence of in- 
dependent random variables (bj)j>o with values in K, then the dynamics of the 
periodic Lorentz gas in the Boltzmann-Grad limit can be described in terms of the 
discrete time Markov process defined as 

(S j+1 ,H j+1 ) = T hj (Hj), j>0. 

Denote Z j+1 = R 2 x S 1 x R+ x [-1, 1] x K n+1 for each j > 0. The asymptotic 
independence above can be formulated as follows: there exists a probability measure 
P on R + x [—1,1] such that, for each j > and ^ <E C{Z n+ \), 
(H) 

lim / ^(x,u},rT r (f,uj),h r (^,uj ),h r ,...,b r n )dxduj 

J rZr xS 1 

ty(x,cu, r, h, b , . . . , h n )dxdujdP (T 7 h)dfi(b ) . . . dfi(h n ) , 



where /i is the measure defined in (2.27). 

This scenario for the limiting dynamics is confirmed by the following 

Theorem 2.7. Let f m be any continuous, compactly supported probability density 
on R 2 x S . Denoting by R[8] the rotation of an angle 9, let F = F(t,x,u>, s,h) be 
the solution of 
(2.30) 

{d t +u-V x - d s )F(t, x, u, s,h) = J 2P(2s, h\h')F(t, x, R[0(h')]u, 0, h!)dti , 



F(0,x,uj,s,h) = f in (x,oj) 



rOC /> 1 

/ / P{ Tl h\h')dh'dT , 
J2s 7-1 



where (x, u,s,h) runs through R 2 x S 1 x R^x] — 1, 1[, and 8(h) = tt — 2 arcsin(/i). 
Then the family (/ r )o<r<i defined in (1-4) satisfies 

(2.31) f r ~* J J F(-,-,-,s,h)dsdh in L°°(R + x R 2 x S 1 ) weak-* 



(IS 



r -> 0+. 



Let us conclude this presentation of our main results with a few remarks. 

Equation (2.30)-(2.31) was proposed first in [7], under some additional decor- 
relation assumption left unverified — specifically, assuming (H). Then, Marklof- 
Strombergsson provided a complete, rigorous derivation of that same equation in 
[19], without any additional assumption, thereby establishing the theorem above. 

The main novelty in this description of the Boltzmann-Grad limit of the periodic 
Lorentz gas is the fact that it involves a Markov process in the extended phase space 
R 2 x S 1 x R + x [—1, 1]. In addition to the space and velocity variables x and w 
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that are usual in the classical kinetic theory of gases, this extended phase space 
involves two extra variables: i.e. s, the (scaled) time to the next collision and h, 
the impact parameter at that next collision, as additional coordinates describing 
the state of the moving point particle. To the best of our knowledge, the idea of 
using this extended phase space (and particularly the additional variables s and h) 
appeared for the first time in our announcements [13, 7]. 

3. Main Results II: Dynamical Properties of the Limiting Equation 

The present section establishes some fundamental mathematical properties of 
equation (2.30). For simplicity, we henceforth restrict our attention to the case 
where the space variable x varies in the flat 2-torus T 2 = R 2 /Z 2 . 

3.1. Equilibrium states. As is well-known, in the kinetic theory of gases, the 
equilibrium states are the uniform Maxwellian distributions. They are characterized 
as the only distribution functions that are independent of the space variable and 
for which the collision integral vanishes identically. 

In equation (2.30), the analogue of the Boltzmann collision integral is the quan- 
tity 

r-l 

2P(2s, h\h')F(t, x, R[6(ti)]w, 0, h')dti + d s F(t, x, w, s, h) . 

-l 

On the other hand, the variables (s, h) play in equation (2.30) the same role as the 
velocity variable in classical kinetic theory. 

Therefore, the equilibrium distributions analogous to Maxwellians in the kinetic 
theory of gases are the nonnegative measurable functions F = F(s, h) such that 

-d s F(s,h)=J 2P(2s,h\h')F(0,h')dh' , s>0,-l<h<l. 

Theorem 3.1. Define 

r-+co r-l 

E(s,h) := / / P(T,h\h')dh'dT. 

J2s J-1 

1) Then 

E(0,h) = l, -Kh<l 

and 

+ 00 r-l />+00 r-l r-l 

/ E(s,h)dhds= / / / \SP(S,h\ti)dhdti dS = 1. 
J-i Jo J-i J-i 

2) Let F = F(s, h) be a bounded, nonnegative measurable function on R + x [—1,1] 
such that s i— > F(s, h) is continuous on R + for a.e. h e [—1, 1] and 

-d s F(s,h)=J 2P(2s,h\h')F(0,h')dh' , s > , -1 < h < 1 , 

lim F(s,h) = 0. 

Then there exists C > such that 

F(s, h) = CE(s, h), s > , -1< h< 1 . 

3) Define 2 

p(t)= lim l{(^ w )e(^n[0,l] 2 )xS 1 |2rr r (x,c)>OI 



r^0+ \(z r n [0,1] 2 ) X s 1 



2 The existence of this limit, and an explicit formula for p are obtained in [2]. 
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Then 

j E(s,h)dh = -2p'(2s), s>0, 

and 

f 1 1 

/ E(s, h)dh ~ ass^ +00 . 

Notice that the class of physically admissible initial data for our limiting equation 
(2.30) consists of densities of the form 

F in (x,co,s,h) = f in (x,uj)E(s,h) 

- see Theorem 2.7. In other words, physically admissible initial data are "local 
equilibrium densities" , i.e. equilibrium densities in (s, h) modulated in the variables 
(x,u>). 

Before going further, we need some basic facts about the evolution semigroup 
defined by the Cauchy problem (2.30). The existence and uniqueness of a solution 
of the Cauchy problem (2.30) presents little difficulty It is written in the form 

F(t,;;-)=K t F in , t>0, 

where (K t )t>o is a strongly continuous linear contraction semigroup on the Banach 
space L 1 (T 2 xS'x R + x [—1, 1]). It satisfies in particular the following properties: 

1) if F in > a.e. on T 2 x S 1 x R + x [-1, 1], then, for each t > 0, one has K t F m > 
a.e. on T 2 x S 1 x R + x [-1, 1]; 

2) for each t > 0, one has K t E — E; 

3) if F m < CE (resp. F m > CE) a.e. on T 2 x S 1 x R+ x [-1, 1] for some constant 
C, then, for each t > 0, one has K t F m < CE (resp. K t F m > CE) a.e. on 
T 2 x S 1 x R+ x [-1,1]; 

4) for each F m <= i x (T 2 x S 1 x R+ x [-1, 1]) and each t > 0, one has 

[[[[ K t F m dxdudsdh = [[[[ F m dxdwdsdh; 

J J J JT 2 xS 1 xR + x [—1,1] J J J JT 2 xS 1 xR + x [—1,1] 

5) if F m e C(T 2 x S 1 x R + x [— 1, 1]) is continuously differentiable with respect to 
x and s, i.e. V x F m and d s F m e C(T 2 x S 1 x R + x [-1, 1]), then, for each t > 0, 
one has K t F m , d t K t F m , V x K t F m and d s K t F m e C(T 2 x S 1 x R + x [-1, 1]), and 
the function (t, x, lj, s, h) i-> K t F m (x 1 co, s, h) is a classical solution of the equation 
(2.30) on R^ x T 2 x S 1 x R; x [-1, 1]. 

All these properties follow from straightforward semigroup arguments once (2.30) 
is established. Otherwise, the semigroup (K t )t>o is constructed together with the 
underlying Markov process in section 6 of [19] — see in particular Propositions 6.2 
and 6.3, formula (6.16) and Theorem 6.4 there. 

3.2. Instability of modulated equilibrium states. A well-known feature of the 
kinetic theory for monatomic gases is that generically, local equilibrium distribution 
functions — i.e. distribution functions that are Maxwellian in the velocity variable 
and whose pressure, bulk velocity and temperature may depend on the time and 
space variables — are solutions of the Boltzmann equation if and only if they are 
uniform equilibrium distribution functions — i.e. independent of the time and 
space variables. In other words, the class of local Maxwellian states is generically 
unstable under the dynamics of the Boltzmann equation. An obvious consequence 
of this observation is that rarefied gas flows are generically too complex to be 
described by only the macroscopic fields used in classical gas dynamics — i.e. by 
local Maxwellian distribution functions parametrized by a pressure, temperature 
and velocity field. 
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Equation (2.30) governing the Boltzmann-Grad limit of the periodic Lorentz gas 
satisfies the following, analogous property. 

Theorem 3.2. Let F be a solution of (2.30) of the form 

F(t,x,w,s,h) = f(t,x,u)E(s,h), with f e C 1 (7 x T 2 x S 1 ) , 

where I is any interval of R + with nonempty interior. Then f is a constant. 

Thus, the complexity of the equation (2.30) posed in the extended phase space 
T 2 x S 1 x R + x [—1, 1] cannot be reduced by postulating that the solution is a local 
equilibrium, whose additional variables s and h can be averaged out. 

As in the case of the classical kinetic theory of gases, this observation is important 
in the discussion of the long time limit of solutions of (2.30). 

3.3. H Theorem and a priori estimates. In this section, we propose a formal 
derivation of a class of a priori estimates that includes an analogue of Boltzmann's 
H Theorem in the kinetic theory of gases. 

Let h be a convex C 1 function defined on R + ; consider the relative entropy 

H h (fE\E) ■= 

(h(/) - h(l) - h'(l)(/ - 1)) (t, x, w, s, h)E(s, h)dxdudsdh . 



/T 2 xS 1 xR+x[-l,l] 

The most classical instance of such a relative entropy corresponds with the choice 
h(z) = z\nz: in that case h(l) = while h'(l) = 1, so that 

H zXnz {fE\E) = [[( (flnf-f+l)(t,x,uj,s,h)E(s,h)dxdujdsdh. 

J J J T 2 xS 1 xR + x [— 1,1] 

Theorem 3.3. Let F = F(t,x,uj, s,h) in C 1 (R+ x T 2 x S 1 x R + x [-1,1]) be such 
that 

< F(t, x, u, s, h) < CE(s, h) , (t, x, oj, s, h) e R+ x T 2 x S 1 x R + x [-1, 1] , 
and 

{d t +u-V x - d s )F(t, x, w, s,h) = J 2P(2s, h\h')F(t, x, R[9{ti)]w, 0, h')dh' , 
with the notations of Theorem 2.7. Then H h (F\E) € C 1 (R+) and 



^-H b (F\E)+ I D h (F/E)(t,x)dx = 0, 
at Jrp2 



/ T 2 

where the entropy dissipation rate Dh is given by the formula 
Dh(f)(t,x) = 

llll 2P(2s, h\ti) (h(f(t, x, R[6(h')]u>, 0, h'))-h(f(t, x, u, s, h)) 

J J J JS 1 xR+ x [— 1,1] x [— 1,1] v 

- h'(f(t, x, oj, s, h) (f(t, x, R[0(ti)]w, 0, h') - f(t, x, oj, s,h))^j dujdsdhdh' . 
Integrating the equality above over [0,i], one has 

H h (F\E)(t)+ f f D h (F/E)(T,x)dxdT = H h (F\E){0) 
Jo Jt 2 

for each t > 0. Since h is convex, one has 

H h (F\E) > and D h (F/E) > 0, 
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and the equality above entails the a priori estimates 
f < H h (F\E)(t) < H h (F\E)(0) , 



I \ D h (F/E)(t,x)dxdt < H h (F\E){0) . 
Jo Jt 2 

That Hh(F\E) is a nonincreasing function of time is a general property of Markov 
processes; see for instance Yosida [25] on p. 392. 

3.4. Long time limit. As an application of the analogue of Boltzmann's H The- 
orem presented in the previous section, we investigate the asymptotic behavior of 
solutions of (2.30) in the limit as t — > +oo. 

Theorem 3.4. Let f in = f in {x,u) G L°°(T 2 x S 1 ) satisfy f in {x,uj) > a.e. in 
(x,u>) e T 2 x S 1 . Let F be the solution of the Cauchy problem (2.30). Then 

F(t,;;;-)^CE 

in i°°(T 2 x S 1 x R + x [-1, 1]) weak-*, with 

C =h\l f n (x,uj)dxduj. 

J JT 2 xS 1 

3.5. Speed of approach to equilibrium. The convergence to equilibrium in the 
long time limit established in the previous section may seem rather unsatisfying. 
Indeed, in most cases, solutions of linear kinetic models converge to equilibrium in 
a strong L 2 topology, and often satisfy some exponential decay estimate. 

While the convergence result in Theorem 3.4 might conceivably be improved, the 
following result rules out the possibility of a return to equilibrium at exponential 
speed in the strong L 2 sense. 

Theorem 3.5. There does not exist any function $ = $(£) satisfying 

= o(t~ 3/2 ) as t -> +oo 

such that, for each f m G i 2 (T 2 x S 1 ), the solution F of the Cauchy problem (2.30) 
satisfies the bound 

\F(t, ') - (/ m )-B||L2( T 2 xS i xR+x [_ 1 !]) 

< $(t)\\F(0, ; ; ; ')||l 2 (T 2 xS 1 X X [— 1, 1] ) 

for each t>0, with the notation 

(^) = h II <p{x,uj)dxduj 

J JT 2 xS 1 

for each<Pe L X (T 2 x S 1 ). 

By the same argument as in the proof of Theorem 3.5, one can establish a similar 
result for initial data in L P (T 2 x S 1 ), with the L 2 norm replaced with the LP norm 
in (3.1), for all p e]l,oo[; in that case $(t) = o(i-( 2p - 1 )/ p ) is excluded. 

The case p — 2 discussed in the theorem excludes the possibility of a spectral 
gap for the generator of the semigroup K t associated with equation (2.30) — that 
is to say, for the unbounded operator A on L 2 (T 2 xS'x R + x [—1, 1]) defined by 

Af(x, co, s, h) = (u-V x - d s )f{x, to, s,h)-j 2P(2s, h\h')f(x, R[6(h')]w, 0, h')dh! , 
with domain 

D{A) = {f G L 2 (T 2 xS 1 xR+x[-l,l])\(oj-V x -d s )f G ^(T^S^R+x [-1, 1])} . 



(3-1) 
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4. An ergodic theorem with continued fractions 



4.1. Continued fractions. Let a e (0, 1) \ Q; its continued fraction expansion is 
denoted 

(4.1) a=[0;a 1 ,a 2 ,as,...} = 



ai + 



1 



1 



d2 



03 + . . . 



Consider the Gauss map 
(4.2) 



T : (0,l)\Q3i4-- 

x 



e (0,1) e (0,1) \Q. 



The positive integers 01,02,03,... are expressed in terms of a as 



(4.3) 



01 



and 



1 



->n-l 



Q: 



n > 1 . 



The action of T on a is most easily read on its continued fraction expansion: 
(4.4) T[0; 01,02,03, . . .] = [0; a 2 ,a 3 ,a 4 , . . .} . 

We further define two sequences of integers (p n ) n >o an d (q n )n>o by the following 
induction procedure: 

Pn+1 = a n Pn + Pn-1 , PO = 1 , Pi = , 

9n+i = a n q n + q n - X , q = , Oi = 1 . 



(4.5) 



The sequence of rationals (^-)„>i converges to a as n — > oo. Rather than the usual 
distance I— — a I, it is more convenient to consider 



(4.6) 

Obviously 



d„ := \q n a-p n \ = (-1)™ (<z„a-p„), n>0. 



(4.7) d n+ i = —a n d n + d n -\ , rfo = 1 , d\ = a . 

We shall use the notation 

o„(a), p„(a), g„(a), d n (a), 

whenever we need to keep track of the dependence of those quantities upon a. For 
each a £ (0, 1) \ Q, one has the relation a n (Ta) = a n+ i(a), which follows from 
(4.4) and implies in turn that ad n (Ta) — d„ + i(a), for each integer n > 0, by (4.7). 
Therefore, 



n-l 



(4.8) 



d n (a) = ]jT k a, n>0. 



While a n (a) and d n (a) are easily expressed in terms of the sequence {T k a)k>a 1 the 
analogous expression for q n (a) is somewhat more involved. With (4.5) and (4.7), 
one proves by induction that 



(4.9) 
Hence 



q n (a)d n+1 (a) + q n+1 (a)d n (a) = 1 , n > . 



i \j i \ 1 d n+ i(a)d n (a) 

q n+1 (a)d n (a) = 1 - ^g„(a)d n _i(a) . 

a„(Q!)a„_i(a) 



so that, by a straightforward induction 

(4.10) q n+1 (a)d n (a) = ^(_ r i ^(»)W n > . 
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Replacing (dk{ot))a<k<n+i by its expression (4.8) leads to an expression of q n+ i(a) 
in terms of T k a for < k < n. 

4.2. The ergodic theorem. We recall that the Borel probability measure dG(x) = 
InSTTx on (0, 1) is invariant under the Gauss map T, and that T is ergodic for the 
measure dG(x) (see for instance [16]), and even strongly mixing (see [21].) 
For each a G (0, 1) \ Q and e G (0, 1], define 

(4.11) N{a,e) = M{n>0\d n {a) < e} . 
Lemma 4.1. For a.e. a G (0, 1), one has 

N(a,e)~^ln±, e -+ 0+ . 

See [6] (where it is stated as Lemma 3.1) for a proof. 
We further define 

(4.12) Sn (a,e) = ^-, n>0, 
for each a G (0, 1) \ Q and e G (0, 1]. 

Theorem 4.2. For m > 0, let f be a bounded measurable function on (R + ) m+1 ; 
then there exists L m (f) G R independent of a such that 

1 F da 

l^YJ^j J f( S N(a,e)(a),dN( a ,e)~l(a),---,$N(a, t )~ m (a))— ~> L m (f) 

and 

f(-l) N ^f(S N(a , e) (a), . . .,S N{a>e) - m (a))^ -+ 
for a.e. a G (0, 1) as rj — »■ + . 

Proof. The proof of the first limit is as in [6], and we just sketch it. Write 

f 1 de 

\ /(5iV(a,e)(a),*iV(a,e)-l(a)! • • • ,^N(c,e)-m(a)) — 

J n e 

N(a,e)-1 „ dn _ l(a ) & 

= V] / /(^iV(a,e)(a),*iV(a,e)-l(a),---)^Ar(a,e)-m(a)) — 

„=1 ^n(a) £ 

JV(Q ' e) - l(Q) , , , „ , , , , ,,de 



/■aw(a,e)-i W 

+ / /(^JV( a ,e)(a),^JV(a,e)-l(a),---,^JV(a,e)-m(a)) 



c 



Whenever d n (a) < e < d n -\(a), one has N(a, e) — n so that, for a.e. a G (0, 1), 

/ /(^iV(a,e)(a))^(a,e)-l(a)) ■ ■ j *iV(a,e)-m(«)) ~ 

iV(o^)-l ^ 

u),..., 5„_ m (a)) — +0(1). 



= /(<5„(a),(5„_i(a) 

n=l J <(«) 



Substituting p — dn ^ in each integral on the right hand side of the identity above, 
one has, for n > m > 1 



/" 



/ f(6„(a),6 n -i(a),..., Vn—m 

(«))- 

•/ <*»(<*) 6 



(a)/ dn _ 1( a) 7 v ' r«-ia ' • • • ' nr=i T»-™+k-laJ p 

= F m _ 1 (T n ~ m a). 
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with the notation 

dp 
P 



Thus, for a.e. a € (0, 1), 

l n (l/ v ) j f( S N(a,e)(^)^N(a,e)-l(^), ' • ' ^ N (a.e)-m(a)) ~ 
1 N(a - V) ^ ( 1 \ 



We deduce from BirkhofT's ergodic theorem and Lemma 4.1 that 

N(a,r))-1 j 

(4A3) £ Fm-ifT^a) ^ L m (f) := ^ j o F m ^(x)dG{x) 

for a.e. a € (0, 1) as r\ — > + , which establishes the first statement in the Theorem. 
The proof of the second statement is fairly similar. We start from the identity 



e 



/ (-l) Ar(a:e) /(<5w(a,e)(a),<5 A r (Q:e )_i(Q;), . . . , £jV(a,e)-m(aO) 

V (-1)"/ /(<y n (a),<5 n _i(a),...,5„_ m (a))-+0(l) 
„=i Jdnia) e 

AT(a,rj)-l 

= £ (-l) n F rn _ 1 (T"- m a) + 0(l). 



Writing 



JV(a,?7)-l 



n—m 



(F m ^(T 2k - m a)-F m ^(T 2k+1 - m a))+0(l) 

m/2<fe<(JV(a,7j)-l)/2 

we deduce from BirkhofT's ergodic theorem (applied to T 2 , which is ergodic since 
T is mixing, instead of T) and Lemma 4.1 that, for a.e. a <G (0, 1) and in the limit 
as i] — » + , 

1 ]T (F m _i(T 2fc - m a)-F ro _i(T 2fc+1 - m a)) 



^ f 1 (F m _ 1 (x)-F m _ 1 (Tx))dG(x)=0 
Jo 



10 

since the measure (iG is invariant under T. 

This entails the second statement in the theorem. □ 

4.3. Application to 3-obstacle configurations. Consider to e S 1 such that 
< UJ2 < wi and W2/001 Q, and let r G (0, 5)- 

The parameters (A(u, r), B(co, r),Q(co, r),a(cu, r)) defining the 3-obstacle config- 
uration associated with the direction to and the obstacle radius r are expressed in 
terms of the continued fraction expansion of W2/0J1 in the following manner. 
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Proposition 4.3. For each ui £ S 1 such that < 0J2 < Wi and ^2/^1 4- Q> an d 
each r <E (0, |) one has 



A(w,r) 
B(u>,r) 



_ dmpht) ( a ) 

e 



dN(a,e)-l(o<) 



djv(a, £ )(a) 



djv(a,€)(a) 



Q(w,r) = eq N{a ^) 



where 



a = 



U>2 
LO-i 



and e = 



2r 
uji 



By the definition ofN(a,e), one has 

< A(co, r) , B(w, r) , Q(co, r) < 1 . 

This is Proposition 2.2 on p. 205 in [6]; see also Blank-Krikorian [1] on p. 726. 
Our main result in this section is 

Theorem 4.4. Let K = [0, l] 3 x {±1}. For each F e C(K), there exists C(F) e R 
such that 

1 f 1/2 dr 
^l^y I F(^(a;,r),B( W ,r),Q( W ,r),(7( W ,r)) T £(F) 

for a.e. uj G S 1 smc/i i/iai < W2 < ^1, in the limit as 77 — > + . 
Proof. First, observe that 

B, g, a) = B, Q) + B, Q) 

with 

F± (A, B, Q) = i B, Q, +1) ± F(A, B, Q, -1)) . 

By Proposition 4.3, one has 

F±(A(u,r),B(uj,r),Q(uj,r)) 

1 - <5jV(a,e)-l(oO 



= F± I 1 - 5jv(a,e)(a), 1 - 5jV(a,e)-l(") 



1 



<^JV(a,€)(a) 
-rfjV(a,e)-l(a)?iV(a,e)(a) 



Sn 



N 



For each m > 0, we define 

B± ^l-5jv,l-5iv-i- 
Observe that 

aTa < \ for each a e (0, 1) \ Q , 
so that, whenever n > m + 1, 



<W(a,e)-l(a) 

fm,±{$N, ■ ■ ■ , SN-m-l) ■ 



J = (iV-m) + 



q n {a)d n -i(a) - ^ (-1) 



n _ •d n (a)d n _i(a) 



j=n—m 



dj(a)dj_i(a) 



< 



d 7t (a)d Tl _i(a) 

77,-777 — 2 (a) 



< ] J T fe a • T fe_1 Q! < 2-" 1 " 1 . 



k—n—m — 1 
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Likewise, whenever n > m + 1 and n > I + 1, 
(4.14) 

^ n _j dn(a)d n -i(a) ^ , iAn _ ? - d n (a)d w _i(a) 



, l J (a) 

1—TI / J \ ' J V / '1 = 1) 



(-If 



dj(a)dj-i(a) 

j—n—m 'j 



< 2 



■ min(/,m) 



Since 5j\r( a e )_i(a) > 1 by definition of AT (a, e), one has 
(4.15) 

|/m,±(<W(a,e)(aO, • • • ,5jV(a,e)- m -l(aO) _ /;,±(^JV(a,e)(a), . . . , £jV(a,e)-!-l(aO)| 

< i o±(2- min ( m '')). 

where p± is a modulus of continuity for F± on the compact [0, l] 3 . 

The inequality (4.14) implies that 
(4.16) 

.1/4 



I /" i/4 dr 
— J^ F ± (A(u,r),B(w,r),Q(uj,r))- 



In(l/4»j) 

1 /-l/2«i ^ e 

~ ln(l/47?) / 2 /m>±(^JV(a,e)(a),---,^JV(a,e)-m-l(Q!))— <P±(2~ m ~) 

upon setting a =002/^1- Moreover, by Theorem 4.2, for each 77* > 

T—f 7T / /m,+ (5jV(a,€)(a))---)^JV(a,e)-m-l(a)) ^Wl(/m, + ) 

(4 17) ln W*/^A 6 

1 /"'* rfr 

y (-lf^/ m ,-(^ M («),-,^ a , f )- m -iW) 7 ^0 

a.e. in a G (0, 1) as 77 — > + , and the inequality (4.15) implies that 

IWi(/ m ,+)-%(/i, + )l<ft(r min(m ' l) ), l,m>0. 

In other words, (£ m +i(/ro,+))m>o is a Cauchy sequence. Therefore, there exists 
LeR such that 

(4.18) £m+i(/m,+) -> i as n -> 00 . 

Putting together (4.16), (4.17) and (4.18), we hrst obtain 

L m +i(fm,+) - 2p(2- m - 1 ) 
1 f 1/4 dr 

n^0+ ln(l/4?7) J v r 



and 



<L m+1 (f m , + ) + 2p(2- m - L ), 

-2p(2- m ~ 1 ) 

1 f 1/4 dr 

* M wm I ^ r)F - r) ' r) ' Q( "' r)) ~ 

1 /-V4 dr 

" i m + RTT^) I ^r)F-(A{w,r),B( U ,r),Q(^r)) T 

< 2p(2- m - 1 ) . 
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Letting m — > oo in the above inequalities, we finally obtain that 

i ^/ fl 1/ V + (^r) )J B( W ,r),0( W ,r))^-,L 

X ' A dr 
<j(u, r)F_ (A(u, r), B(w, r), Q(w, r)) ► 



ln(l/4»7) 7, 

a.e. in w € S 1 with < UJ2 < oj\ as 77 — > + . □ 

Amplification of Theorem 4.4. The proof given above shows that 

£(F)=£(F + ) 

for each F e C(K). 

5. Computation of the asymptotic distribution 
of 3-obstacle configurations: 
a proof of Proposition 2.6 

Having established the existence of the limit £(F) in Theorem 4.4, we seek an 
explicit formula for it. 

It would be most impractical to first compute L m+ i(/ m + ) — with the notation 
of the proof of Theorem 4.4 — by its definition in formula (4.13), and then to pass 
to the limit as m — > 00. 

We shall instead use a different method based on Farey fractions and the asymp- 
totic theory of Kloosterman's sums as in [2] . 

5.1. 3-obstacle configurations and Farey fractions. For each integer Q > 1, 
consider the set of Farey fractions of order < Q: 

F Q :={Z\l<p<q<Q, g.c.d.(p,g) = l}. 
If 7 = 2 < 7' = £7 are two consecutive elements of Fq, then 

q + q' > Q, and p'q — pq' = 1 . 
For each interval I C [0, 1], we denote 

Fq{I)=IC\Fq. 

The following lemma provides a (partial) dictionary between Farey and continued 
fractions. 

Lemma 5.1. For each < e < 1, set Q = [1/e]. Let < a < 1 be irrational, 
and let 7 = | and 7' = £7 6e ffce two consecutive Farey fractions in Fq such that 
7 < a < 7'. TTien 
ft f < a < ^ iften 

?iV(a,e) ( a ) = Q and d N{a.e) {(*) = qOl - p ] 

(11) if^ < a < 2±£, tfien 

to(a, E )(a) = min(o,g/) 

aV(a,£) (a) = qa-p if q<q , and d N(ate) (a) = p' - q a ifq'<q; 
(Hi) if ^-t! < a < y> then 

«jv(a,€)(a) = 9' div(a, e )(a) = p' - q'a . 
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Sketch of the proof. According to Dirichlet's lemma, for each integer Q > 1, there 
exists an integer q such that 1 < q < Q and dist(ga, Z) < Qqrp If ^77 G Tq is 
different from - and ^, then q cannot be equal to q" . For if a < ^ < £77, then 
p" — q"a > -^r(p"q' — p'q") > ^77 > q^j- Thus q is one of the two integers q and 
q'. In the case (i) p' — q'a > e > Qzp[ so that q 7^ q', hence q = q. Likewise, in the 
case (hi) qa — p > e > so that q = q' . In the case (ii), one has 

0< ^(1-qe) = ^j{qp' - pq 1 - qe) < qa - p < e 

since q <Q<\. Likewise < p' — q'a < e, so that q is the smaller of q and q' . □ 

In fact, the parameters {A{uj, r), B(w, r), Q(co, r)) can be computed in terms of 
Farey fractions, by a slight amplification of the proposition above. We recall that, 
for each w € S 1 such that < cj 2 < u>i and cj 2 /wi is irrational, one has 

Q(w,r) — eqN( a e)i a ) 7 w ith a — — and e = — . 
Under the same conditions on uj, we define 

JJluj, r) := , again with a = — and e = — , 

e oji lji 

and Q'(uj,r) in the following manner. 

Let Q = [1/e] with e = and let 7 = - and 7' = 2. be the two consecutive 
Farey fractions in Tq such that 7 < a < 7'. Then 

(i) if 2 < a < we set 

Q'(uj,r) := eg' ; 

(ii) if pLz! < a < 2+1 we set 

W q' — — q 1 

Q'(u),r) := em&x(q,q') ; 

(hi) if < a < 2l , we set 

Q'(u,r) := tq. 

With these definitions, the 3-obstacle configuration parameters are easily expressed 
in terms of Farey fractions, as follows. 

Proposition 5.2. Let < r < j and weS 1 be such that < w 2 < CJi and oj 2 /u;i 
is irrational. Then 

A(u,r) = 1 - D(w,r) , 
6(w,r) 



B(w, r) = 6(cj, r) 

with 

b(u,r) = 



D{u,r). 



Q{uj,r) - l + Q / (w,r)D(w,r) 



Q(w,r) 

Sketch of the proof. We follow the discussion in Propositions 1 and 2 of [2] . Con- 
sider the case | < a < ^r^, and set e? = <ja — p > and d' = p' — q'a > 0. 
According to Proposition I of [2], eB(uj 7 r) = e — (pk — <ZfcCti), with the notation 
Pk = p' + fcp ct qt = q' + kq for fc G N* chosen so that 

p' + kp — e p' + (k — l)p — e 

q' + kq a q' + (k - l)q ' 
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In other words, k is chosen so that d' — kd < e < d! — (k — l)d, where d = qa — p 
and d' = p' — q'a. That is to say, 



'e-d 1 ' 




'l-d'/e 


d 




D 



and 



-k 



n , d 1 -kd , d' d , 

B = 1 = 1 \-k- = b- 

e e e 



D 



with 6=1-^. Since qp'-pq' = 1, one has d' = (l-q'd)/q, so that 6 = Q ~ 1 q ' £ ' . 



The other cases are treated similarly. 



□ 



5.2. Asymptotic distribution of (Q, Q' , _D). As a first step in computing C(F), 
we establish the following 



Lemma 5.3. Let f e C([0, l] 3 ) and J = [a_, a+] C (0, 1). Then 



f{Q{u, |ewi), Q'(u, \ewi), D(u, iewi))l W2 / WieJ 



[04] 



f(Q,Q',D)d\(Q,Q',D) 



as e — > + 7 where A is i/ie probability measure on [0, l] 3 given 6y 



12 



(lo<Q : Q'<llQ+Q'>ll 0<I)< ig3 + lQ<Q'lig3 < _ D< i) 



dX(Q,Q',D) = 
dQdQ'dD 



Q 



The proof of this lemma is based on the arguments involving Kloosterman's sums 
to be found in [2]. 

Proof. For < e < 1, set Q = [1/e]. Observe that 



Q(uj, |ewi) - — Q(w, |ewi) 



< 1 Q{u,\eu{) 1 



Q+l eQ "Q + l' 



< 1 g^, < 1 



Q+l eQ 



Q+l ' 



since |ewi) < eQ and Q'(w, |ewi) < eQ while 

|D(w, |ewi) - eQD(uj, ±ewi)| < e - - Q 



< e < 



Q 



since D(cu, \eoJ\) < 1. 

Since / is continuous on the compact set [0, l] 3 , it is uniformly continuous; let p 
be a modulus of continuity for /. Then 



/ f(Q(u, |ewi), Q'(cj, |ewi), D(u, |ewi)) — 
Jsi_ " " w l 



(5.1) 



< 3p(l/Q) 



as e 



0+. 
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Henceforth, we seek to prove that 



(5.2) 



J gi /^Q(w,5ewi),^Q'(w,3ewi),eeD(w,|ewi)^ 

-+\J\I f(Q,Q',D)dX(Q,Q',D) 
J[o,i] 3 



as e — > + . Without loss of generality, by an obvious density argument we restrict 
our attention to the case where the test function / is of the form 

f(Q,Q',D)=g(Q,Q')h(D), 

with g e C([0, l] 2 ) and h G C([0, 1]). Then 

= E // ff(!,£W(9«-p))<fa 

p/<i^ a {J) J q w ^ 

+ E / P ^(|4)MQ( P '-^)Ma 



(5.3) 



p/qeF Q (J) q 



+ E fp'-e 1 9<«'-9 ( 3 ) h{Q{qa - p))da 



p/qEFaiJ) q' 

+\\9h-\\i 







_ P 


+ sup 


p' 


_ P 


sup 












9' 


3 


p'/q'er Q {J) 


q' 


q 



= I + II + III + IV + 0(1/Q) 



where - < ^ are consecutive elements of To . Notice that the O term above is the 

q q 

contribution of the cndpoints of J not being elements of Tq . 
Define 



H(z) = / h(y)dy , for each zeR. 
Jo 



Thus 



p/qeF Q (J) 
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One has 



Hi Q 



Hi Q 



,i-«/e 



< 



1 



likewise 



while 



Hence 



with 



Hi Q 



1 - eq' 



Hi Q 



< 



q'Q(Q + l) ' 
1 



Q 



q Q(Q+1) ' 



\H(1) - H(eQ)\ < \\h\\ L oo\l - eQ\ = \\h\\ L ~e 



e 



< \\h\\ L ~e. 



p/ q er Q (J) 

|*/|<IMU~|W| L ~ E 



p/qer Q (J) 

Applying Lemmas 1 and 2 in [2] shows that 



q'Q(Q + l) 



E 9 

p/q£F a (J) 



1 l\ 1 H ( l - q / Q \ 

Q' QJ Qq \ q'/Q J 



dQdQ' 
>i — q — as Q -> oo . 



while 



I J\ ■ £ jj g{Q, Q')H {^) 



\RA < \\g\\L~\\h\\ L ~± ■ |J| ■ 4, w>1 ^^ = o(i/Q) 



so that 



(5.4) / | J| • 4, ^ S (Q, Q')# (^) 1 Q+Q '>i^ as Q -+ oo . 
By a similar argument, one shows that 

H ~* ' J| ^ ^ // i » 5(Q '' Q) " lQ+Q,>1 ^ ' 

(5.5) IH^\J\-£ ll o ^9(Q,Q')(H(l)-H(^))l Q<Ql l Q+QI>1 ^, 

IV^\J\-js ff 9(Q'> Q) ( H (!) - # (*#)) lQ'<QlQ+Q'>i^ 



as Q — > oo. 

Substituting the limits (5.4)-(5.5) in (5.3) shows that 



dQdQ'dD 



g{Q,Qi)h(D)\ Q+QI>1 \ 0<D< ^^ 

(0,1) 3 



+kl ■ ^ Ijj 9(Q, Q')h(D)l Q<Q ,l Q+Q 
as e — > 0, which establishes (5.2). 



dQdQ'dD 

'>lJ-i_Q< D <l Q 
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On account of (5.1), this proves the convergence announced in Lemma 5.3. □ 



5.3. Asymptotic distribution of (A, B,Q). Next we compute the image the 
probability measure A under the map (Q, Q', D) M> (A, B, Q) defined in Proposition 
5.2. In other words: 

Lemma 5.4. Let g e C([0, l] 3 ). Then 

-+ / g(A,B,Q)dv(A,B,Q) 
J[o,i] 3 

as e — > + 7 w/iere ^ is f/ie probability measure on [0, l] 3 given 6?/ 

dAdSdQ 



12 



<MA 5 ^ Q) = — 1o<a<i1o<b<i-a1 <q< 



1 - A 



Proof. We first compute the image $*A of the probability measure A in Lemma 5.3 
under the map 

$ : [0,1] 3 ^ [0,1] x R x [0,1] 

(Q, Q', D) i ^ (A, 6, Q) = (1 - £), ^i±S^, Q) 
The Jacobian of $ is 



detV$(Q,g',Z?) = 



n ^Q.'-p i 



-1 



21 o 
Q u 



D I- A 



0' 



A straightforward computation shows that 



This expression can be put in the form 



l-^<6<0 + 1 fc>olb>2-A-i 



d$*A(A,6,Q) - 
dAdbdQ 



\ dAdbdQ 
J 1-A ' 



dAdbdQ 



(5.6) <2<I>*A(A, &, Q) = ^ 3 -lo<A,Q<l(lfceA 1 (A,Q) + l6eA 2 (A,Q)) 

where 

A 1 (A,Q) = (A-^,(1-^)A0) 
A 2 (A,Q) = ((2-A-i)V0,(l- A)) 

The probability measure ^ is the image of ^A under the map 
* : [0, 1] x R x [0, 1] — > [0, l] 3 

(A, 6, Q) i ^ (A, B, Q) = (A, b - [^(1 - A), Q) 

In other words, 

i / a T-) /^\\ 12 -. -. dAdbdQ 

df{A,±S,Q) = ^-lo<A,Q<llo<B<l-A 1 _^ 



(5.7) 



• £ (1a i( a,q)(S + n(l - A)) + 1a 2 ( A ,q)(B + n(l - A))) 



_ 12 



^j-M (A, if, Q)i-0<A,Q<li-0<B<l-A—^-^- 
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where 



(5.8) 



M(A, B,Q) = #{n eZ\B + n(l - A) e A^A, Q)} , 
+ #{n g Z | B + n(l - A) e A 2 (A, Q)} . 



Whenever u, w </ Z, one has 



#{™ eZ|me( u ,D)} = ([«] - [ u ])+ . 



Hence, for a.e. A,Q e [0, 1] 



#{n E Z\B + n(l - A) E A X (A, Q)} 



(1- A/Q) AO- B 
1-A 



A-A/Q-B 



I- A 



since 



(1- A/Q) AO- B A-A/Q-B (1 - A) A (A/Q - A) 



I- A 



I- A 



I- A 



> 0. 



For a.e. A, Q E [0,1], if Q < A 



(1- A/Q) AO -B 



I- A 

l-A/Q-B 



A- A/Q - B 



I- A 



I- A 
A-A/Q-B 



A-A/Q-B 
1-A 



I- A 
A-A/Q-B 
I- A 



Likewise, for a.e. A,Q E [0, 1] 



#{n E Z | B + n{\ - A) E A 2 (A, Q)} 



lq>a 



l-A/Q-B 



I- A 



(2- A- 1/g) VO-5 



1-A 



observe that, whenever Q < A, one has 2 — A — 1/Q < 2 — A — 1/A < 0, so that 



l-A/Q-B (2- A- 1/g) V0-B 1 - A/Q 



1-A 



1-A 



1-A 



<0, 



while, if g > A, 



1 - A/Q - B (2 - A - 1/g) V0-B 



1-A 1-A 

(i - A/g) a (A -i)(i- i/g) 

1-A 



> 0. 
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Therefore 



#{n e Z I B + n(l - A) e Ai(A, Q)} 
e Z | B + n(l - A) e A 2 (A, Q)} 



(1- A/Q) AO -B 



1-A 



A- A/Q-B' 
1- A 



+ 



1 - A/Q - B 



1-A 
= 1q<a + 1q>a 



(2- A- 1/Q) VP -B 

(1- A/Q) AO- B 
1-A 

(2- A - 1/Q) VO-5 



= 1 + 1q>a| 





' -B " 






1- A_ 





1-A 

(2- A - 1/Q) VO-B 



Now, if Q < 



(2 - A - 1/Q) V - B 



1-A 



2-A 





' -B ' 




1-A 


that 





M{A, B,Q) = #{n e Z | B + n(l - A) e Ai(A, Q)} 
+ #{n eZ|B + n(l-4)e A 2 (A, Q)} 

"2-A-l/Q-B 





' -B " 











1 + 



If < Q < 1 and < B < 1 - A, one has 



1-A 



so that 



-(1 - A) < 2 - A - 1/Q - B < (1 - A) 
-B 



1-A 
2-A-l/Q-B 

i - A 
2- A- 1/Q -B 
1-A 



= -1, 

= 0, if 2-A- £> 1/Q, 
= -1, if2-A-B<l/Q. 



Therefore 



(5.9) 



M(A,B,Q) = 1 + 1 Q>3 ^ 



B 



1-A 



2-A- 1/Q- B 
1-A 



~ 1 ~~ 1 Q>^ 1 2-A-B>l/Q - 1q <5 _i_^ 

whenever B/0, and this establishes the formula for v in the lemma. 



□ 



5.4. Computation of C{F). With the help of Lemma 5.4, we finally compute 
C(F). 

Proof of Proposition 2.6. Let F e C(K), and set 

F+(A, B,Q) = \ (B(A, B, Q, +1) + B(A, B, Q, -1)) ; 
obviously B+ e C([0,1] 3 ). 
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By Lemma 5.4 



J[0,1\ 



F+(A,B,Q)dv(A,B,Q) 



as 77 — > + , for each 77* > 0. 

On the other hand, by Fubini's theorem 



Mv*/V) J v 



J' (J F+(i4(w,|a«;i),S(w,ia«;i),Q(a;,|ea;i)) 



efcj \ de 



/ f 1 r wa F 

is; \In(j?*/?7) J Wir) /2 



L(.1U'.,).^.,).^.,-))-)^. 

r / u>f 



By Theorem 4.4, the inner integral on the right hand side of the equality above 
converges a.e. in lj £ S\ to C(F+). By dominated convergence, we therefore 
conclude that 

C{F+)= [ F + (A,B,Q)dv(A,B,Q). 
Finally, according to the amplification of Theorem 4.4 

C{F)= f F + {A,B,Q)dv{A,B,Q) 

J[o,i] 3 

or, in other words, 

C(F)= I F(A,B,Q,a)dv{A,B,Q,a) 

J[0,l] 3 x{±l} 

where 

d/i(A, B, Q, cr) = \dv{A, B, Q) ® (5 CT=+1 + S a= -{) . 

□ 



6. The transition probability: a proof of Theorem 2.1 

6.1. Computation of P(S, h\h'). We first establish that the image of the proba- 
bility v defined in Lemma 5.4 under the map 

[0, l] 3 3 (A, B, Q) ^ T AtBiQi+1 (h') € R+ x [-1, 1] 

is of the form P+(S, h\h')dSdh, with a probability density P + which we compute 
in the present section. Let / G C C (R + x [—1,1]); the identity 



L 



f(S, h)P+(S, h\h')dSdh = f f(T AtBiQi+1 (ti))dv(A, B, Q) 

R+x[-l,l] J[0,1] 3 
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defining P+(S, h\h') is recast as 

f f(S,h)P+{S,h\h')dSdh 
Jr + x[-i,i] 

= / f(Q,ti -2(l-A))l 1 _ 2A<h , <1 du(A,B,Q) 

+ / / ( ^ Q i-A B \ h ' + 2(1 - B)) l- 1<hl< - 1+2B diy(A, B, Q) 
J[o,i] 3 v ' 

+ I f { ^ Q i B I A \ h' + 2(A - B)) l-i+2B<h><i-2Adv(A, B, Q) 
J\o,i] 3 v ' 



I + II + 11 1 . 



In I, we first integrate in B, since the argument of / does not involve B. Ob- 
serving that V (2 — A — q) < B < I — A, we have to distinguish two cases, namely 
< Q < Tf^j an d < Q < 1; which leads to 

I=^J [f(Q, h'-2(l-A))l {1 _ h , )/2<A<1 (l 0<Q< ^ + z ^l^ <Q<1 ) dQdA . 

In the right-hand side of the equality above, we set S = Q and h = h' — 2(1 — A), 
to find 

II h)1 - l<h<h,<1 ( 1 o<s< 1+(fc /_ h)/2 + |(^yl T ^<s<i) d Sdh . 
With the notation r) = ^\h — h'\, we see that I can be put in the form 



I = ^JI f(S, h)l- 1<h<hl<1 (s v l 0<s<rh + (1 - S)1__l 



m75 <5<1 7 Sri 



dSdh 



f(S, /i)l_i< fc <v<ilo<s<i(S»?) A (1 - 5)^ 
As for II, given A, we first make the substitution 

{S , h )= (^q^l,h' + 2(1- B)) 

whose Jacobian is 



d(S, h) 

Hence 



l-A u 
l-A Z 



2(1 -B) _h-b! 
l-A ~ l-A 



, 1 dSdh \ i A 

A+il-H'y* < S <lht L A<^<1 \h-h'\ 



The inner integral is recast as 



/ 1 i iSi i 1 .. h-h' ^-.dA 

J l-A+(h-h>)/2 <a< l-A A<^— <i 

= ls>ol-l<h'</Kl / ^-1-1/ S<A<\-1/ S+(h-h') /2^- A< h^h!_dA 

Jo 

= ls>ol-K/i'</i<lls<l(l + f) + 



+ ls>0l-l<h'</Klls>l(^ L 2* L + i - 1) + 
= ls>ol-K/i'<h<l(^- - |4 - 1|)+ , 
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so that, eventually, we find 

= jj f(S, h)l s>0 l-i< h <<h<i(Sv - |1 - 

Now for III. Here we make the substitution 

(S,h) = C-^ A - B \ h' + 2(A-Bj) 

whose Jacobian is 



d(S, ft) 



A-B 
' l-A 

Q o 

l-A Z 



2(A-B) _ h-h' 
l-A ~ l-A 



This Jacobian vanishes for A = B; therefore, we further decompose 

in = iih + iih 

with the notations 



Hh = I f f 1 "^"^ ,^ + 2(^4 - S)) 1_ 1+2B<V<1 _ 2A 1 A<B ^(A, B, Q) , 
J[o,i] 3 v 7 

///2 = /" / A ) )/t / + 2(yl _ B) \ i_ 1+2B</i/<1 _ 2A i a>b ^ (j 4 5 Bj Q) . 

J[0,1] 3 v 7 

We begin with 77 ii . With the substitution above, one has 

IHl = J iyj J /(^' ^) 1 2A-K(/i-ft')/2<Al/i'<l-2Al-l + 2A<h 

li/(i-A)<s<i/(i-A+(ft-/ l ')/4) 1 - Kh</i'<irsz77f^ <k4 

= f(S,h)^J ^-2A-l<(h-h') /2<A^-h' <1-2A^- -1+2 A<h 

^-l/(l-A)<S<l/(l-A+(h-h')/A)dA^ l-Kh</i'<l jfz^J • 

The inner integral is 

/ l2A-l<(/i-h')/2<Al/i'<l-2Al-l+2A<hll/(l-A)<S<l/(l-A+(h-/i')/4)^ 

Jo 



A) 1 ^<2 1 + i( ft - h ') 1 ' 4< 3( 1 -' l ') 1 ' 4< 2"( 1+ ' l ) ll_ l + 3(' l_ ' l ') <A<1_ | rfj4 

= ((I + ±(ft - ft')) A |(1 - h') A |(1 + ft) A (1 - i) 

-(1- 1 +I(ft-ft'))V0) + 

= ((l-^ + /i'|-i(ft'-/i))A(l-i)-(l-i-i(/ i '-/ l )|)VO) + . 

Alternatively, using the relations (a + c)V(fe+c) = aVfe+c, (a+c)A(6+c) = aAfe+c, 
(—a) A (—6) = — (a V &) and aV6 + aA& = a + fe, we find that 

((i -\\h + h'\- \{h' - ft)) A (1 - i) - (1 - 1 - ±(ft' - ft)) V 0)+ 
= (H - II* + a (|(ft' - h) - §) + l-±(ft' -h)-{l- \{h> - ft)) V I + i)+ 
= ((1 - ±(ft' - ft)) A i - (i + ±|ft + ft'|) V (1 - \{h' h))) + 
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SO that 



Hh = $ J J f(S,h)^J l2A-K(h-h')/2<Alh><l-2Al-l+2A<h 

X ll/(l-A)<S<l/(l-A+(h-h>)/4)dA S j l-l<h<h'<l (IfJh' l 

= % J J f(S, h)l- 1<h<hl<1 ((S - \S(ti — ft)) A 1 

- { l-l S (h'-h))V(§ + \S\h + h'\)] dS ' lh 



' + S\h-h>\ ' 

With the notation £ = ||ft + ft'|, we recast this last expression as 

dSdh 



I"i = £ f(S,h)l_ 1<h<h><1 ((S-lSr 1 )Al-(l-lS V )V^S+lSO) + 



The computation of III2 is fairly similar. Starting with the same substitution 
as for Ilh, we obtain 

J iyj J f(^^)^-2A-l<(h-h')/2<A^-h'<l-2A^--\+2A<h 
ll/(l-A+(/ t -fe')/4)<a<l/(l-A)l-Kfe'</Kl [ff-h'^ ^A 

fff^'ty^f 1 2A-K(h-h>)/2<Alh><l-2Al-l+2A<h 

j a\ , dSdh 
1 l/(l-A+(h-h')/4)<S<l/(l-A) dA i—Kh'<h<l 1 



|ft-ft'| 



The inner integral is recast as 



L 







l2A-K(/i-h')/2<Alh'<l-2/ll-l+2A<ftll/(l-A+(/i-h')/4)<S<l/(l-A)^ 



1 

Q 1 ^(h-h^KA 1 A<^ + j(h~h' ^ AK^il-h') 1 AK^l + h) 1 1-^ <A<l-^ + j(h-h' ) dj,L 
1 

1 i(/ l -ft')<^ 1 ^<5-il'»+'»'l+5('i-'»') 1 i-^<^<i-|+3('J-'i') d " 4 
= ((1 _ 1 + + l {h - h >)) A (i_ § + i(ft- h'))- (l-ijv^- ft')) + . 

As in the case of , we have 

((! -\\h+h'\ + \{h h 1 )) a (1 - 1 + \{h - h 1 )) - (1 - 1) v - //))+ 
= (1 + - ft') + (-1 _ + /,'d A (-i) - (1 - 1) v - h')) + 

= {i + \{h-h')-(\ + \\h + h'\)y\ 

_(l--i)-i(ft-ft') + (l-l)A±(ft-ft')) + 
= ( | - | (ft - ft') - (| + I|ft + h'\) V (I) + (1 - I) A ±(ft - ft'))+ 
= ((1 - ±(ft - ft')) A (i + \{h - ft')) - (I + ±|ft + ft'|) V ±) + 
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so that 

Ilh 



— -p- J {^j j f(S, h)l 2 A-l<(h-h')/2<Alh><l-2Al-l+2A<h 



1 | dSdh 

i l/(l-A+(/ l -/ l ')/4)<S<l/(l-A)-l— Kh'<h<l |7^77| 



\dA 



12 



In other words 



f(S, h)l- 1<h , <h<1 ((S - \S(h - h 1 )) A (1 + \S(h - h 1 )) 

dS 

s\h- 



- (iff+iS | fc + fc1) vi) 



m 2 =£ JJ f(s,h)i- 1<h<h , <1 ((s-±s v )A(i^s v )-(±s^sovi) + ^ . 



P+(S, h\ti) = ^ r l_ 1<ft</l , <1 lo<s<i(5? ? ) A (1 - 5) 



+ Tcv: 1 - K?i'<'i<ils , >o(5'?7 - I 1 - 



+ l^l-Kk<v<ils>o ((S-±S77)Al-(l-±Sr,)V(±S+±SC))_ 



Summing up the contributions /, II, and 77/2, we find that 
6 

_6_ 

_6_ 

+ ^l-i<fc'<fc<ils>o ((S-|S»7)A(l + iS» 7 )-(|5+i5C)Vl) + . 

Observe that 

l s> o(^ - |1 - = l^^^Sr, - |1 - 5|) . 

On the other hand, if 

(S-^)A1-(1-^)V(|5+|5C)>0 

then (S'-iS'??) - (I-^-Stj) > so that S > 1; likewise ±-(±£+±50 > so that 
S < 2. Hence 

ls>o((S-|S»7)Al-(l-i5f ? )V(|S+i50) + 
= 1 1<5<2 ((5-IS^)Al-(l-|5» 7 )V(|S+i50) + • 
By the same token 

{S-±Sr])A(l + ±Sri)-{±S+±SC)Vl > 
implies that (5— \Sr}) — 1 > 0, so that S > 1, and hence 

l S >o((5-|S»7)A(l + iS >7 )-(iS+iSOvl) + 
= l s >i ((S-|5»j)A(l + |S»7)-(|5+|SC)Vl) + • 

Finally 

6 

+ ^ 1 - 1</l ' </l<ll T^<5<T^^ - I 1 - 5 D 

+ -^;l-i<h<fc'<ili<s<2 ((S-±S V )Al-(l-±Sr 1 )V(±S+±SO) + 



P+(S, h\ti) = -^;l_i< h <v<ilo<s<i(S»7) A (1 - 5) 



+ ^l-i<v</,<il5>i ((5-i^)A(l + i^)-(l5+i5C)Vl) + . 
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By formula 2.24, 
so that 

/i|/0 = | h\ti) + P+(S, -h\ - h')) , 

thereby leading to formula (2.11). 

6.2. Proof of the simplified formula (2.13). Assume that h > \h'\ so that 

T)=\{h-h') ( h=(( + V ) 

C = \(h + ti) \ h' = {c-v) 

and denote 

L = S-±Sr], M = ±S+±S(, N=±Sr). 
Whenever S > 1, the reader will easily check that (2.12) can be written as 

^NP(S,h\h') = {1 + N-L)+ + (LA(1 + N)-MV1)+ + (LA1-MV(1-N)) + . 

Since L> M, the expression above can be reduced after checking the three cases 
L < 1, l<L<l + A r andi>l + iV. One finds that 

(1 + N - L)+ + (L A (1 + N) - M V 1)+ = (1 + N - M V (L A 1))+ . 

Then 

ZfNP(S, h\ti) = (1 + N - L)+ + {L A (1 + N) - M V 1)+ 
+ (LA1-MV(1-JV))+ 

= {1 + N -MV (LA 1))+ + (LM-MV(1- iV))+ 

if M > 1 + N , 

1 + N -M ifl-N<M<l + N, 
2N if M < 1 - N . 

Since 

M + N =\S{l + h) 

M-N =\S(l + ti) 
the formulas above can be recast as 

' if \S{1 + ti) > 1, 



2tt 2 



NP{S,h\ti) = < 



l-|5(l + /i') if + ti) < 1< ±S(1 + /i) , 



\S{h-ti) if 1 < + 

which holds whenever S > 1. 

On the other hand, if S < 1, the last two terms on the right hand side of (2.12) 
vanish identically so that 

^NP(S, h\ti) = (27V) A (1 - S) + {2N - (1 - S))+ = 2N = \S{h - ti). 

Putting together these last two formulas leads to (2.13). 
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6.3. Proof of Corollary 2.2. As for the statements in Corollary 2.2, observe that 
the symmetries (2.14) and the fact that P is piecewise continuous on R+ x [—1, 1] x 
[1, 1] follow from (2.12). The first identity in (2.15) follows from the fact that, by 
definition, (s,h) i->- P(s, h\h!) is a probability density, while the second identity 
there follows from the first and the symmetry P(s, h\h!) = P(s, h'\h). 

If S > 4 and h, h! €] — 1, 1[ satisfy \h'\ < h, the simplified formula (2.13) implies 
that 

P(S,h\ti) < -j^— ^l 1+fcf< a . 

On the other hand, 1 + h' < § and S > 4 imply that hi < -\, so that h > \h'\ > |; 
therefore h — h' > 1 and the inequality above entails (2.16). 

Starting from (2.16), observe that, because of the symmetries (2.14), 

f f P{s,h\ti)dhdti = 4 ff P(s,h\h')dhdh' 

J-1 J-1 JJo<\h'\<h<l 

24 f f 48 
3c // l-i+h^zdhdti = -— , 

K & JJ0<\h'\<h<l ^ TT -3 



< 



which is (2.17). 



7. Proof of Theorem 3.1 

The existence of the integral defining E follows from the positivity of P(s, h\h') 
and the second identity in (2.15). 

That E(0, h) = 1 for \h\ < 1 follows from the formula defining E and the second 
identity in (2.15): this establishes the first part of statement 1). 

The definition of E and the first formula in statement 1) show that 

-d s E(s,h) = J 2P(2s,h\h')dh' = J 2P(2s, h\h')E(0, h')dh' , s > , \h\ < 1 . 

Using again the positivity of P(s, h\h') and the second identity in (2.15) shows that, 
for each h £ [— 1, 1], the function (s,h') ^ P(s,h\h') belongs to L 1 (R+ x [—1, 1]). 
Therefore E(s,h) — > for each h e [—1,1] as s — > +oo. Thus each function 
F = F(s, h) of the form F(s, h) = CE(s, h) satisfies both conditions in statement 
2). 

Conversely, let F = F(s, h) satisfy the conditions in statement 2), and let := 
F(0,h) for a.e. h e [—1,1]. Integrating both sides of the differential equation 
satisfied by F in s £ R + yields 

(7.1) $(/i) = J Il(h\h')$(h')dti 
since F(s, h) vanishes as s — > +oo, with 

/•oo 

IL(h\h') = / P(S,h\ti)dS. 
Jo 

Multiplying each side of the identity (7.1) by $(/i), and integrating in h £ [—1, 1], 
we see that 

(7.2) J <$>{K) 2 dh = J J Il(h\ti)${h)${ti)dhdti . 

Observe that, by (2.14) 

n(h'\h) = U(h\h') , for a.e. h, ti £ [-1, 1] , 
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and 



J U(h\h')dh' = J Il(h'\h)dh' = 1 , for a.c. h e [-1, 1] 



in view of (2.15). 
Therefore 



J <S>(hfdh = J (^J Ii{h\ti)dh^j<5>{h) 2 dh = J Tl(h\ti)dhJ §{tifdh , 
so that (7.2) becomes 

= J §{hfdh- J j n(h\ti)<P(h)<P(ti)dhdti 

= J J n(h\ti)±($(hf + $(tif)dhdti - J J n(h\ti)<t>(h)$(ti)dhdti 



/ n(h\ti)\(<b(h)-<s>{ti)) 2 dh. 
i J-i 



Since H(h\h') > a.e. in h, h! € [—1, 1] as can be seen from the explicit formula 
(2.18), this last equality implies that 

<P(h) =C a.e. in fee [-1,1], 

for some nonnegative constant C . Therefore 

-d a F(s,h)=C [ 2P(2s,h\h')dh! , and lim f(s» = a.e. in h e [-1, 1] , 

J —I s— ^+oo 



so that 



(-00 rl 

F(s,h) = C / P(r,h\h')dh'dr = CE(s,h), 

J2s J-l 

which proves the uniqueness part of statement 2). 

Now for statement 3); by definition of P(s, h\h'), for each t > 

pi poo pi pi 



f E(s,h)dh= f f f P{t,h\ti)dtdhdti 

J-l J2s J-l J-l 

1 f 1/4 f 1 dr 
= £ 1 ™ f 77^ J f y_ i M2s,+^M-ui](T r (h , ,uj))dh / — 

1 /- 1/4 f 1 dr 
= £ Um J J i l [2s , +co )(2rr r (r r (fe', W ))^' — 

i f 1/4 

= 2 lim - — - / v r ({(x,w) e T+/Z 2 I 2rr r (x,w) > 2s}- 
e^o+ | lne| J e 



" 1/4 ,„ , ^,„o,„ ^ 

r 

where v r is the probability measure on that is proportional to u> ■ n x dxduj. 

Using formula (1.3) in [2], which is a straightforward consequence of variant of 
Santalo's formula established in Lemma 3 of [9], we conclude that 

ri 

E(s, h)dh = -2p'(2s) , s>0, 

-l 

which is the first formula in statement 3). The second formula there is a consequence 
of the expression of p"(s) as a power series in 1/s given in formula (1.5) of [2]. 
Finally, we establish the second formula in statement 1). Indeed 

r-oo pi 



/: 



/ / E(s, h)dhds = / -2p'(2s)ds = p(0) = 1 . 
Jo J-i Jo 
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while the first equality in that formula follows from the identity defining E and 
Fubini's theorem. 

8. Proof of Theorem 3.3 

If F is a solution of (2.30), we set / = F/E, so that, observing that E(0, h') = 1 
for ti € [-1,1], 

E(s, h)(d t + W'V X - d s )f(t, x, uj, s, h) - f(t, x, u, s, h)d s E(s, h) 



Since 



= J 2P(2s,h\ti)f(t,x,R[6(h%0,ti)dti 
d s E(s,h) = - j 2P{2s,h\h')dh' , 



the equation above can be put in the form 

E(s,h)(d t + uj ■ V x - d s )f(t,x,u>,s,h) 

= J 2P(2s,h\h')(f(t,x,R[0(h')],0,ti)- f(t,x,uj,s,h))dti . 

Next we multiply both sides of the equation above by h'(/) — h'(l): the left 
hand side becomes 

(h'(/(t, x, u, s, h)) - h'(l))E(s, h)(d t + uj-V x - d s )f(t, x, u, s, h) 

= (d t +oj ■ V X -8 S ) (E(s, h)(h(f(t, x, w, s, h)) - h(l) - h'(l)(/(i, x, uj, s, h) - 1))) 

+d s E(s, h)(h(f(t, x, uj, s, h)) - h(l) - h'(l)(/(i, x, to, s, h) - 1)) 

= {d t +uj ■ V X -8 S ) (E(s, h)(h(f(t, x, uj, s, h)) - h(l) - h'(l)(/(i, x, uj, s, h) - 1))) 

- J 2P(2s, h\h')(h(f(t, x, uj, s, h)) - h(l) - h'(l)(/(t, x, uj, s, h) - l))dh' . 

Integrating in (ui,s,h) transforms this expression into 

f fOO P 1 

d t / {h(f) - h(l) - ti {l){f -l)){t,x,w,s,h)E(s,h)dhdsduj 

Js 1 Jo J-l 

+ div x / / w(h(f) -h(l) -h'(l)(/- l))(t,x,w,s,h)E(s,h)dhdsduj 
Js 1 Jo J-l 

+ I [ (h(/) - h(l) - h'(l)(/ - 1))) (t, x, uj, 0, h)dhduj 

- I I II 2P(2s,h\h'){h(f)-h(l)-h'(l)(f-l)){t,x,w,s,h)dh'dhdsduj. 
Js 1 Jo 7-i j-i 



Using the relation (2.15) simplifies this term into 

roo i-l 



d t I I I (h(/) -h(l) -h'(l)(/ -l)){t,x,uj,s,h)E{s,h)dhdsduj 
Js 1 Jo 7-i 

p pOO p 1 

+ div x / / / uj(h(f) -h(l) -ti(l)(f -l))(t,x,uj,s,h)E{s,h)dhdsduj 
Js 1 Jo 7-i 

+ / / (h(f)-ti(l)f)(t,x,uj,0,h')dh'duj 
Js 1 J-l 

III 2P(2s,h\h')(h(f) -h'(l)f)(t,x,uj,s,h)dtidhdsduj . 

s 1 Jo 7-i 7-i 
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On the other hand, multiplying the right hand side by h'(f(t, x, ui, s, h)) — h'(l) 
and integrating in (co,s,h) leads to 

p pOO p1 p1 

/ / / / 2P{28,h\h')(f(t,x,R[0{h')]w,O,h')-f{t,x,w,8,h)) 
Js 1 Jo J-l J-l 

x (h'(f(t, x, w, s, h)) - h! {l))dti dhdsdw 

p poo pl p1 

= 2P(2s,h\ti)(f(t,x,R[6(ti)]w,0,ti) - f(t,x,u,a,h)) 

Js 1 Jo J-l J-l 

x h'(f(t, x, uj, s, h))dh' 'dhdsdu 

P pOO p1 p1 

+ h'(l) / / / / 2P(2s,h\h')f{t,x,u>,s,h)dh'dhds(ko 
Js 1 Jo J-i J-i 

-h'(l) / / /(i,x,w,0,/i')^' dw - 
Js 1 J-i 

after substituting ui for R[6(h')]u> in the last integral above. 

Putting together the left- and right-hand sides, we arrive at the equality 

r />oo f-1 

d t / (h(f) -h(l) -h'(l)(f -l))(t,x,w,s,h)E(s,h)dhdsduj 

Js 1 Jo J-l 

+ div x [ [ [ w(h(/) -h(l) -h'(l)(/- l))(t,x,w,s,h)E(s,h)dhdsduj 
Js 1 Jo J-l 

+ f [ (h(/) - h'(l)/)(t, x, oj, 0, h')dh'duj 
Js 1 J-i 

p pCO p1 p1 

- / / 2P(2s,h\h')(h(f) - h'(l) f)(t,x,oj,s,h)dh'dhdsdco 

Js 1 Jo J-i J-i 

p pOO p1 p1 

= 2P(2s,h\ti)(f(t,x,R[6(ti)}uj,0,ti) - f(t,x,w,s,h)) 

Js 1 Jo J-l J-l 

x h'(f(t, x, uj, s, h))dh' dhdsdu 

+ h'(l) [ [ [ [ 2P(2s,h\h')f(t,x,co,s,h)dh'dhdsduj 
Js 1 Jo J-i J-i 

-h'(l) J J f(t,x,u>,0,h')dh'doj. 

All the terms with a factor h'(l) in the integral part of the equality above compen- 
sate, so that the equality above reduces to 

r /*oo />1 

d t / (h(f) -h(l) -ti(l)(f -l))(t,x,w,s,h)E(s,h)dhdsduj 

Js 1 Jo J-l 

p /*oo r-1 

div x / / u){h(f)-h(l)-ti(l){f-l))(t,x,u),s,h)E(s,h)dhdsdu) 
Js 1 Jo J-l 



+ 



1 

h(f)(t,x,ui,0,h')dh'dui 
s J-i 

OO f-1 f-1 

/ / 2P(2s,h\h')h(f)(t,x,ui,s,h)dh'dhdsdco 
s 1 Jo J-l J-l 

OO f-1 f-1 



f- f-QO f-1 f-1 

- 2P(2s,h\h')(f(t,x,R[e(h')]u;,0,h')-f(t,x,uj,s,h) 
Js 1 Jo J-l J-l 

x h'(f(t, x,uj 7 s, h))dh! ' dhdsduo = . 
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Using again the relation (2.15) and the substitution uj n- R[8(h')]uj, one has 



/ / h(f)(t,x,w,0,h')dh'dw 
Js 1 J-i 



rOO pi pi 

/ / 2P(2s,h\h')h(f)(t 7 x 7 uj 7 0,h')dh'dhdsduj 
's 1 Jo J-i J-i 

p pCO pi pi 

= / / 2P(2s,h\ti)h(f)(t 7 x } R[9(ti)}u,0,ti)dtidhdsduj 

Js 1 Jo J-i J-i 

so that the previous identity can be put in the form 

p pOO p 1 

d t / (h(f) -h(l) -ti(l)(f -l))(t,x,u,s,h)E(s,h)dhdsdoj 

Js 1 Jo J-l 

+ div x / / w(h(/) -h(l) -h'(l)(/ -l))(t,x, w, s,h)E(s,h)dhdsdco 
Js 1 Jo J-l 

+ J 1^ J J 2P(2s, h\ti) (h(/) (t, x, R[6{ti)]u>, 0, h') - h(/) (t, x, oj, s, h) 

- (f(t,x,R[0(h')]cj,0,ti) - f(t,x,w,s,h))h'(f(t,x,w,s,h)jdh'dhdsdw = 0. 

Integrating both sides of this identity in x e T 2 , we finally obtain 
d f 

-H h (fE\E)+ D h (f)(t,x)dx = 0, 
dt J T 2 

with Hh(F\E) and -Dh(/) defined as in the statement of Theorem 3.3. 

9. Proof of Theorem 3.2 

Inserting F(t,x,u>, s,h) — f(t,x,u>)E(s,h) in (2.30), one finds that 
E(s,h)(d t +uj- V x )f(t,x,cj) - f(t,x,u)d a E(s,h) 



= J 2P(2s,h\ti)f(t,x,R[9(ti)]u)dti 



so that 

E{s,h){d t +wV x )f(t,x,uj) 



since 



= J 2P(2s,h\h')(f(t,x,R[9(h')}cj) - f{t,x,w))dti 

-d s E(s,h)= J 2P(2s,h\h')dh' , s>0,|ft|<l 
- see Theorem 3.1. 

Integrating in h e [—1,1] both sides of the penultimate equality, we obtain 

J E(s,h)dh(d t +UJ-\7 x )f(t,x,uj) 
r-l r l 

2P(2s, h\ti) (f(t, x, R[0(h')]u) - f(t, x, w)) dti . 



l J-i 



Since 



f 1 1 

/ E(s, h)dh ~ _ 9 9 ass^ +oo 



7T 2 S 2 



by Theorem 3.1, while 

rl ,1 C , 



f f P{S,h\h')dhdh' < — 
J -i J -i (1 
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by (2.17), one has 



L 



E{s,h)dh(d t +Lj-V x )f(t,x,Lu) 



:(d t +UJ ■ V x )f(t,x,co) , and 



1 ,-\ 



2P(2s, h\ti) (f(t, x, R[9(ti)]uj) - f{t, x, w)) dh'dh 



< 



4C" 



(l + 2s)3 
(dt+u-V x )f(t,x,u)=Q, and 

d s E(s, ti) (f(t, x, R[6(h')]u) - f(t, x, u)) dti = 



Hence 



Integrating the second equality in s > 0, and observing that -E| s _ = 1 while 
E(s, ti) — > uniformly in ft G [—1, 1] as s — > +oo, we see that 

d s E(s, ti) (f(t, x, R[9(ti)]oj) - f(t, x, u)) dti] ds 



d s E(s, ti)dsj {fit, x, R[9{h')}oj) - f(t, x, co)) ds 

= J ^ (f(t, x, R[6(ti)]uj) - f(t, x, w)) dti = 

or, in other words 

f(t, x,u>) = ±J f(t, x, R[9(ti)]uj)dti =: x) . 

On the other hand, the first equation implies that 

{d t +uj-V x )4>{t,x) =0 

so that 

a t 0-a Xj = o, j = 1,2. 

Hence is a constant, which implies in turn that f(t,x,ui) = F ^^'^ is a 



constant. 



10. Proof of Theorem 3.4 



Let f in e £°°(T 2 x S 1 ) be such that / in (a;,w) > a.e. in (z,w) € T 2 x S 1 , 
and let F be the solution of the Cauchy problem (2.30). Define Fo(t,x,co,s,h) = 
f m (x,uj)E(s, ft) and set KjF = Fj for j € N, where we recall that (K t ) t > is the 
evolution semigroup associated to the Cauchy problem (2.30). 
Step 1: 

Assume first that V™/ m e L°°(T 2 x S 1 ) for all m > 0, so that F solves (2.30) 
in the classical sense. Then, with h(z) = \z 2 , one has 

n 

H{F n+1 \E)+Y,mFi\E) - H{F 3+1 \E)) = H(F \E) 
j=o 

for each n > 0, and since 

H{F 3 \E) - H(F j+1 \E) = H(Fj\E) - H^F^E) > 
by Theorem 3.3, one has 

H(Fj\E) - H(F j+1 \E) ^ as j +oo . 
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Now, 



H(Fj\E) - H(F j+1 \E) = f f D{K t F 3 / E)dxdt 

Jo Jt 2 

P(2s, h\h')$j(t, x, lj, s, h, ti) 2 dhdtidsdujdxdt 



10 Jt 2 Js 1 Jo JJ[-1A} 2 
with the notation 

$, (t, x,u,s, h, ti) = KtF ^ h) S,k) - KtFj (x,R[e(ti)]u,0,ti) 

for a.e. (t, x, u, s, h, ti) G [0, 1] x T 2 x S 1 x R + x [-1, l] 2 . 

In view of properties 1) and 3) of the evolution semigroup (K t )t>o recalled in 
section 3, one has 

K t Fj(x,u,s,h) in 
°" E&h) ~ llf HW* S1 ) 

for a.e. (x, to, s, ti) € T 2 x S 1 x R + x [—1, 1] and each t > 0, so that, up to extraction 
of a subsequence, 

F jk I E—^F/ E in L°°(T 2 x S 1 x R + x [-1, 1]) weak-* 

as jk — > +oo, by the Banach-Alaoglu theorem. Besides, the estimate (2.16) implies 
that 

AC" 

0<(d t +u-\7 x - d s )K t F ]k < __||f«|| LOO(T2xS1) 

so that, by the usual trace theorem for the advection operator d t + uo ■ \7 X — d s (see 
for instance [8]) 

K t F jk (x, R[6(ti)]u, 0, h')^K t F{x, R[6(ti)]u>, 0, ti) 
in L°°([0, 1] x T 2 x S 1 x [-1, 1]) weak-*. In particular 

$ ifc ^$ in L°°([0, 1] x T 2 x S 1 x R+ x [-1, l] 2 ) weak-* 

with 

$(t, x, u, a, h, ti) = K * F ( x > u >B> h ) _ KtF{Xi R[9 (ti)]uj, 0, ti) 
E(s,n) 

a.e. in [0, 1] x T 2 x S 1 x R + x [—1, l] 2 . By convexity and weak limit 

f [ [ [ [ P{2s,h\ti)$(t,x,u,s,h,ti) 2 dtidhdsdu;dxdt 
o Jt 2 Js 1 Jo J-i J-i 

= [ [ D(K t F/E)dxdt < lim f [ D(K t F jk / E)dxdt = , 
Jo Jt 2 fe ^°° Jo jt 2 

$(t, x, u, s, h, ti) = K t F (x>u>s>h) _ ^ R[9{h , )]oj ^ Qj h>) = 
E(s, n) 

a.e. on [0, 1] x T 2 x S 1 x R + x [— 1, l] 2 . Averaging both sides of this identity in 
ti e [—1,1] shows that 

KtF E^h) H) = * j\KtF(x,R[6(h')]«,,0,h')dh' =: f(t,x,u>) 

for a.e. [0,r] x T 2 x S 1 x R+ x [-1,1], i.e. 

K t F(x, uj, s, ti) = f(t, x, lo)E(s, ti) . 



so that 
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By Theorem 3.2, one has f(t,x,uj) = C a.e. in (t, x,u) £ [0, 1] x T 2 x S 1 for some 
constant C > 0, so that 

F(x, oj, s, h) = CE(s, h) a.e. in (x, w, s, h) £ T 2 x S 1 x R + x [-1,1]. 

Let us identify the constant C. Property 4) of the semigroup (K t )t>o recalled 
in section 3 implies that 

// / / Fj(x,u>, s,h)dhdsdxdu> 

J JT 2 xS 1 JO J-1 

= ff f n (x,uj)dxduj f f E(s,h)dhds= ff f m (x,uj)dxduj. 

JJ^xS 1 JO J-1 ii^xS 1 

Since F Jk / E—^F/ E in L°°(T 2 x S 1 x R+ x [-1, 1]) weak-* as j k -> +oo, 

// f m (x, u>)dxdu> = II 11 F(x,ui, s, h)dhdsdxdu> 

JJ^xs 1 JJt 2 xs 1 Jo J-i 

= 11 II CE ( S > h)dhdsdxdio = 2irC . 

JJt 2 xS 1 Jo J-1 



Hence 



Fjk /E^C=± f in (x,oj)dxdoj 

J JT 2 xS 1 



in L°°(T 2 x S 1 x R+ x [-1, 1]) weak-*. 

Since the sequence (Fj/E)j> is relatively compact in L°°(T 2 x S 1 x R + x [— 1, 1]) 
weak-* and by the uniqueness of the limit point as j — > +oo, we conclude that 



Fj/E^C=± f in (x,uj)dxdoj 

J JT 2 xS! 

in L°°(T 2 x S 1 x R + x [—1, 1]) weak-* as j — » +oo. Thus, we have proved that 
KjFo^CE in i°°(T 2 x S 1 x R+ x [-1, 1]) weak-* 

where 

C =hll f n (x,uj)dxduj 

J JT 2 xS 1 

P P POO P 1 

= / / / / Fo(x,u!,s,h)dhdsdxdu) , 
7T JJt 2 xs 1 Jo J-i 

whenever V™/" 1 £ L°°(T 2 x S 1 ) for each m > 0, with f in > a.e. on T 2 x S 1 . 
Step2: 

The same holds true if < f in £ L°°(T 2 x S 1 ) without assuming that V™.f m £ 
i°°(T 2 x S 1 ) for each m > 1, by regularizing the initial data in the x- variable. 

Indeed, if (Cc) e >o is a regularizing sequence in T 2 such that Q{—z) = ( e (z) for 
each z £ T 2 , one has 

P P POO P 1 

II II {K t {(e* x Fo) — K t Fa)(x,uj,s,h)(f>(x,uj,s,h)dhdsdxdGJ 

JJt 2 xS 1 Jo J-1 

= 11 I I K t (Fo)(x,oj,s,h)(£ e * x <f>— (f>)(x,u>, s,h)dhdsdxcLo 

JJt 2 xS 1 Jo J-1 

because K t commutes with translations in the variable x. Since 

\\K t F /E\\ L oo( T 2 xS i xR+x [_ 1A }) = ||i r o/-E||L~(T 2 xS 1 xR + x[-l,l]) 

for alH > 0, and 

Ce*x </>-<f>^ Oin L X (T 2 x S 1 x R + x [-1,1]) , 
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we conclude that 

K t (Ce * x F ) - K t F ^0 in L°°(T 2 x S 1 x R + x [-1, 1]) weak-* 
as e — > + uniformly in t > 0. Since we have established in Step 1 that 

Kj{Ce*xF )^ If f n {x,u)dxdLO 
J JT 2 xS 1 

in L°°(T 2 x S 1 x R + x [— 1, 1]) weak-* for each e > 0, we conclude by a classical 
double limit argument that 

Kj{F )^CE with C= ± 1 1 f n (x,u})dxdu} 

JJ^xS 1 

in L°°(T 2 x S 1 x R+ x [-1, 1]) weak-* as j -> +oo. 
Summarizing, we have proved that 

K i F ^^{^JJ J J F (x,uj,s,h)dhdsdxdu?j E 

in L°°(T 2 x S 1 x R + x [-1, 1]) weak-* as j -> +oo. Replacing F with K t F for 
each t € [0, 1] and noticing that 

fi 

-FT t F (x, w, s, h)dhdsdxduj 



T 2 xS x Jo 



= // / / Fq(x,u>, s,h)dhdsdxdu> 
JJt 2 xS 1 Jo J-i 



we conclude that 



K t F ^±-(H I I F (x,u,s,h)dhdsdxdw) E 

\J J T 2 x S 1 JO J-l / 

in L°°(T 2 x S 1 x R+ x [-1, 1]) weak-* as t -> +oo. 

11. Proof of Theorem 3.5 

The argument used in the proof is reminiscent of the one used in [12, 14]. 
Assume the existence of a profile <&(i) such that the estimate (3.1) holds. There- 
fore, for each initial data f m E L 2 (T 2 x S 1 ), the solution 

F(t,;;;-)=K t (f in (;-)E(;-)) 

of the Cauchy problem satisfies 

\\F(t, •, •, •, •)IU 2 (T 2 xSixR + x[-14]) < ||(.r™)^||L 2 (T 2 xSixR + x[-l,l]) 

+ $(t)||F(0, •, •, •, •)IU 2 (T 2 xSixR+x[-14]) 

= (/ m )v / 27r||£ ! ||L 2 (R + x[-i,i]) 

+ <i) ( i )ll/ m |U 2 (T 2 xS 1 )||^||L 2 (R + x[-l,l]) • 

Assume f in > a.e. on T 2 x S 1 ; then F > a.e. on R + x T 2 x S 1 x R + x [-1, 1], 
so that the right-hand side of the equation satisfied by F is a.e. nonnegative on 
R+ x T 2 x S 1 x R + x [-1, 1]. Thus F > G a.e. on R + x T 2 x S 1 x R+ x [-1, 1], 
where G is the solution of the Cauchy problem 

f {d t +u-V x -d s )G{t,x,uj,s,h) = 0, 

\ G(0, x, u, s, h) = / m (x, w)E(s, h) . 

The Cauchy problem above can be solved by the method of characteristics, which 
leads to the explicit formula 

G(t, x, w, s, h) = f m (x - tco, uj)E(s + t,h), 
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for a.e. (t,x,u,s,h) £ R+ x T 2 x S 1 x R+ x [-1,1]. Thus 

'i ') ') ■)IU 2 (T 2 xS 1 xR+x[-l,l]) 

> II^C*) ') '7 ') OIU^T^SixR+xI-l.l]) 

,oc ,1 \ 1/2 



= ||f n |U 2 (T 2 xSi) 1 1 E(a, hfdhd^j 

> 11/^11^(^x8!)^ (jH (J^E(s,h)dh) ds^j 



1/2 



Therefore, if there exist a profile satisfying the estimate in the statement of 
the theorem, one has 



1 /2 

'f (/>»*)*) £ ^'^"H,* 



for each t > 0, and for each f in £ L 2 (T 2 x S 1 ) s.t. f in > a.e. on T 2 x S 1 . 
Let p £ C(R 2 ) such that 

P>0, supp(p) C (-i,i) 2 , 

and set 

Pe(x) : =P(^) 

for each e £ (0, 1). Define p e as the periodicized bump function 

Pe{x) := Pe{x + k). 
kez 2 



Clearly 



|Pe||z,!(T 2 ) — l|Pe||i 1 (R. 2 ) — II P\ L 1 (R 2 ) 



llPe || i2 (T 2) = ||p e ||L 2 (R 2 ) = £||p||l 2 (R 2 ) • 

Choosing f m (x,u) :— p € (x) in the inequality above leads to 

2 \ 1/2 

„ I IIPIUMR. 2 ) , *,^\ fKu ™n 

,1]) 



2 1/2 

J^E(s,h)dh) dA < (^ pgg + $(*)) V2PIU 2 (R + x[-i 



and, letting e — > + , we conclude that 

2 \ 1/2 



(/ (/ 1 jE(s ' /l)d/l ) ds J " ^WH £ IU 2 (R + x[-i,i])- 

That = o(t~ 3//2 ) as £ — > +oo is in contradiction with statement 3) in Theorem 
3.1, which implies that 

2 \ I/ 2 

E(s,h)dh] ds) 



V37T 2 t 3 / 2 



as £ — > +oo, while statement 1) in the same theorem implies that 

l|£|li 2 (R + x[-i,i]) < J J ^E(s 7 h)E(0,h)dhds = J E(s,h)dhds = l. 
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